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Abstract 


Analysis of meteorological, chemical and microphysical data from the airborne SUCCESS mission is reported. Careful 
analysis of the complex DC-8 flight pattern of May 2, 1996 reveals 19 linear or nearly linear flight segments within six main 
geographical areas, which we have analyzed. Significant mountain wave activity is revealed in the data from the MMS and 
MTP instruments on the DC-8, which resembles previous observations of mountain wave structures near Boulder, Colorado. 
Strong mountain-wave-induced upwelling downwind of the Rockies is noted. Turbulence is also noted in regions of the 
mountain wave consistent with overturning near the tropopause. Zonal winds recorded on the ER-2 are shown to be consistent 
with mountain wave breaking at or near critical levels in the stratosphere, consistent with the strong turbulence reported by the 
pilot during the ER-2 flight. These observations have been supported with spectral analyses and modeling studies. “Postcasts” 
of mountain wave activity on May 2, 1996 using the Naval Research Laboratory Mountain Wave Forecast Model 
(NRL/MWFM) predicts both strong mountain wave activity near the tropopause (as measured by the DC-8) and strong 
mountain-wave-induced turbulence in the stratosphere (as encountered by the ER-2). Two-dimensional simulations of fluid 
flow over topography reveal similar isentropic structures to observations. 


Mountain wave-induced upwelling downstream of the Rockies generated wave clouds that were sampled during this flight. 
Analysis of in situ data from the DC-8 shows that the clouds were thin and formed by mountain wave-induced adiabatic 
cooling. A formal parcel advection methodology has been developed and applied to accurately model the production and 
distortion of vertically thin cloud and contrail layers by such gravity waves. The methodology has been used to extend the 
NRL/MWFM model to forecast mountain wave-induced cooling. The new model has been used to simulate stratospheric ice 
cloud formation by mountain waves on synoptic scales, as well as the effect of cooling by random spectra of gravity waves on 
equilibrium aerosol microphysics. The enhanced NRL/MWFM forecasts will play a major role in future NASA aircraft 
campaigns, such as the SAGE III Ozone Loss and Validation Experiment (SOLVE). The parcel-based theory has also resulted 
in unanticipated spin-off research, such as insights into the formation of gravity wave spectra in the atmosphere. 
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1 . Work and Results 


1.1 Mountain Wave Encounter of May 2, 1996 

1.1.1 Mountain Wave Structure from MMS and MTP Instruments on the DC-8 

On May 2, 1996, the DC-8 flew a complex flight pattern aimed at sampling air within a mountain wave cloud over 
the Rocky Mountains near Boulder, Colorado. Strong velocity oscillations were recorded on the aircraft as it intercepted 
mountain waves. Plate l 1 shows a three-dimensional perspective of the DC-8 flight path on this day (shown in red). 
Superimposed in white is the raw 1Hz vertical velocity activity measured by the new DC-8 Meteorological Measurement 
System (DC-8 MMS) [Chan et at., 1998]. Ground-level color surface plots show topographic elevation. It is immediately 
evident that enhanced small-scale velocity activity was measured by the DC-8 MMS over the Rockies. It is also clear that the 
complex flight geometry of the DC-8 on this day complicates any analysis of the sampled wave activity (a full 3D animation 
of this plot is available on the web 2 ). 

To tackle this, a computer program was developed to read in arbitrary data from the instruments on the various 
aircraft platforms deployed during SUCCESS, and then interactively identify “regular” (i.e., nearly linear) flight segments that 
might be amenable to an analysis of wave properties. Results for the May 2 flight are shown in Figure 1. Black curves show 
(from top to bottom) time series of aircraft longitude, latitude, altitude and true air speed during the flight. Blue curves, 
labeled 1-19, show segments of the flight that were approximately linear in longitude, latitude and altitude. We chose not to 
include aircraft air speed in the classification since, as evident from the figure, it remains approximately constant once the 
DC-8 has reached cruise altitude. For test purposes, we also chose one flight segment (#6) with a small curved portion to its 
latitude-longitude trajectory, to assess the effect of curved airplane motion on derived parameters. This curve is marked in red 
in Figures 1 and 2. 

Figure 2 shows three different planar sections of the three-dimensional flight geometry shown in Plate 1. Flight 
segments 1-19 from Figure 1 are superimposed on the latitude-longitude section (top plot) of Figure 2, along with colored 
shading depicting the underling topographic elevations. This grouping illustrates that flight segments 1-19 cluster within 6 
mam geographical “flight legs,” labeled sequentially A-F and shown in orange in Figure 2. Details of these flight segments are 
given below in Table 1 . 


Flight 

Leg 

Flight 

Segment #’s 

Longitude-Latitude 

Transects 

Cruise 

Altitude 

Geographical Location 

Direction 
(modulo 180) 

i 

1 

40-4 IN, 1 00-1 05 W 

10.6 km 

Kansas 

NW,WNW 

B 

2,3 

-41.5N, 1 06-1 08 W 

11.8 km 

Southern Wyoming 

WNW 

C 

4,5 

40.2-4 1.5N, -106W 

11.8 km 

Wyoming-Colorado 

NNW 


6 


11.2 km 



iilllllilil 

7 


10 km 

Wyoming-Colorado, slightly 



8,11,12,16 

39.6-4 1.3N, -105W 

11.8 km 

downwind (east) of Rockies 

NNW 


15 


12.4 km 




9,10 

39.8N, 105-107.5W 

11.8 km 

Across Rockies, approximately 
west of Boulder, Colorado 

WNW 


13,14 


11.2 km 


17 


11.8 km 



in ijiiii 

18 

39-40N, 98-104.5W 

11.2 km 

Colorado-Kansas 

WNW 


19 


7.3 km 




Table 1: Dominant flight legs (labeled A-F) during the DC-8 flight of May 2, 1996 (see Figures 1 and 2). 


Legs A and F were essentially cruising patterns to/from the airfield from/to the Rocky Mountain region. Leg B was a 
sortie into Wyoming looking for wave clouds. Legs C and D were transects approximately parallel to the ridge axes of the 
dominant topography. Conversely, leg E was oriented roughly perpendicular to the dominant mountain ridge axis (see Plate 1 
and/or Figure 2) which, based on standard mountain wave generation ideas [e.g., Bacmeister et ai, 1994; Leutbecher and 
Volkert, 1996], should efficiently profile any larger-scale mountain waves that were forced by flow over the larger-scale 


full page figures such as these are collated in Appendix D, to aid uninterrupted flow of the text, 
point your browser to the following URL: http://uap-www.nrl.navy.mil/dynamics/html/2may96.html 
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quasi-two-dimensional topography here. However, smaller-scale quasi-three-dimensional topographic features may force a 
spectrum of waves at various orientation angles, which would be profiled effectively during any of the flight legs B-E. 

Figure 3 shows isentropes derived from Microwave Temperature Profiler (MTP) data during the four flight 
segments of flight leg E across the mountain ridge (shown in red). The cruise altitude of each flight segment is also shown 
with a colored curve on each plot. We note isentropic structure consistent with flow over this underlying topography, and 
subsequent penetration of forced mountain waves into the upper troposphere and stratosphere. 

Figure 4 shows the three-component velocities, temperature and turbulent energy dissipation rates (TEDRs) 
measured in situ by the MMS during each of these flight legs (with 3 point smoothing applied). TEDR values were 
computed from the high frequency 1-5 Hz MMS data using the method described by Chan et al. [1998]. We note that the 
flight segment pairs at the same altitude (segments 9 and 10, segments 13 and 14) yielded very similar measurements, 
indicating a quasi-stationary mountain wave feature. We also note that the vertical velocity data (Figure 4c) are “noisier” 
and dominated by shorter scales than the horizontal velocity and temperature data (Figures 4a, 4b, and 4d, respectively). 

The TEDR data are noiser again, being dominated by highly intermittent short-scale turbulence events. 

To study the large-scale mountain wave structures more effectively, these data were smoothed using a 65 point 
(~15 km) running average. These smoothed profiles are shown in Figure 5, which “cleans up” the vertical velocity and 
TEDR data in Figures 5c and 5e, respectively. In the zonal velocity and temperature traces (panels a and d, respectively), we 
see clear sinusoidal peaks and troughs with a horizontal wavelength -0.5° (around 43 km), in reasonable agreement with 
earlier estimates [e.g., Dean-Day et al, 1998], In the meridional winds (panel b) and temperatures (panel d), we also see 
evidence of a tonger scale wave with a wavelength -2° (-170 km). We also note order-of-magnitude increases in TEDRs at 
around 105.5°W, where large dips in the isentropes are observed in Figure 3. Similar enhancements in TEDRs are also seen 
at -107°W. 

1.1.2 Wave Cloud and Microphysics Observations from the DC-8 and ER-2 

Wave clouds were also encountered on this flight, principally along flight leg D. Figure 6 shows total aerosol 
concentrations, volumes and surface areas measured by the MASP 3 instrument on this day [Baumgardner and Gandrud, 
1998]. Figure 7 shows ice water content measured by the CVI 4 instrument (top panel) on the DC-8 [Gerber et al., 1998]! 
These data reveal large aerosol volumes and water/ice occurrences on flight segments 8,11,12 and 16, which were all part of 

■ flight leg D (see Figure 1) and all occurred at a 

height of -11.8 km (see Table 1). Note that during 
flight segments 6, 7 and 15, also part of flight leg D, 
these instruments did not record significant aerosol 
or water/ice content. The cruise altitudes of these 
latter flight segments were 10, 11.2 and 12.4 km, 
respectively (see Table 1). This indicates a thin 
vertical cloud layer located at -11.8 km, inferences 
supported by the cloud lidar returns taken from the 
ER-2 during May 2, 1996 (see Figure 8 opposite). 
From Figure 3, this cloud altitude of -11.8 km lies 
just a few hundred meters below the tropopause. 

Flight leg D lies approximately north-south 
at 105°W, somewhat downstream (eastward) of the 
dominant topography (see Figure 2). From Figure 3, 
this region is seen to be at the isentropic ascent phase 
of the mountain wave near the tropopause, yielding 
Figure 8: Backscattered power from the downward-scanning CLS adiabatic cooling and thus increased likelihood of 
(cloud hdar system) on the ER-2, May 2, 1996. cloud particle condensation at -11.8 km. Similar 

processes have been simulated using NRL/MWFM 
simulations of mountain waves in the Arctic stratosphere [Carslaw et al., 1998a, 1998b, 1999], where parcel theory is used 
to predict mesoscale cooling due to wave-induced upwelling [Eckermann and Bacmeister, 1998; see also section 1.4], 
Microphysical simulations of the wave clouds measured here support the notion that wave-induced upwelling and adiabatic 
cooling triggered ice cloud formation via heterogeneous freezing along flight leg D [e.g., Jensen et al., 1998; Weinheimer et 
al., 1998] (see also section 1.4). 


3 multiangle aerosol spectrometer probe 

4 counterflow virtual impactor 
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1.2NRL/MWFM “Postcasts” for SUCCESS Mission Days 

Figure 9 plots maps of “postcast” mountain wave activity for the Boulder region on May 2, 1996. These maps were 
generated using the Naval Research Laboratory’s Mountain Wave Forecast Model (NRL/MWFM 5 ) interfaced to National 
Center for Environmental Prediction (NCEP) assimilated winds and temperatures for the day (daily global NCEP data on a 
5°x2° grid). 

Figure 9c shows postcast peak mountain wave displacements at 200mb, which is within the DC-8 cruise altitude 
range of ~1 1-12 km. The DC-8 flight track on May 2, 1996 is also shown. We see that NRL/MWFM predicts large mountain 
wave displacement amplitudes near the eastward end of flight leg E, as observed in Figures 3 and 4. Figure 9d shows postcast 
normalized momentum flux divergence between levels, a measure of wave dissipation and local turbulence production (0 = 
no dissipation, 1— maximum wave dissipation). We see that NRL/MWFM predicts some turbulence on the westward end of 
flight leg E. Some evidence of such turbulence exists in TEDRs inferred from the MMS in Figure 5e. 

Recent extensions of the NRL/MWFM model have been made so that mountain-wave-induced mesoscale cooling 
can be forecast by the model [see Eckermann and Bacmeister, 1998; see section 1.4], since it is this wave-induced effect that 
is of most relevance to microphysical model predictions. The peak vertical displacement amplitudes of -1-2 km in Figure 9c 
translate into peak mesoscale coolings of -10 K in mountain- wave updrafts near the tropopause, which are sufficient to drop 
air below the frost point and trigger the ice clouds observed on flight leg D [Jensen et al., 1998], 

The ER-2 flew a similar flight pattern through the stratosphere on this day, occasionally overflying the DC-8 in 
efforts to observe from above any subsequent contrailing from the aircraft or effects of aircraft exhaust on the mountain wave 
cloud (the ER-2 flight track is plotted in Figures 9a and 9b). This ER-2 flight was “highly turbulent” according to pilot notes 6 . 
Figure 9b shows normalized momentum flux divergence at lOOmb (-16 km) simulated by NRL/MWFM. We note large 
predicted mountain-wave-induced turbulence production at this level that frequently intercepts the ER-2 flight path, consistent 
with the turbulence reported by the pilot during this flight. The reasons for this intense turbulence production can be noted 
from the navigational zonal winds recorded on the ER-2 during this flight, shown in the top panel of Figure 10. We note zonal 
wind values near zero and even some westward values, as opposed to the strong eastward winds in the troposphere. Such 
vanishing or reversing zonal winds in the presence of zonally-aligned mountain waves are characteristic signatures of the 
wave overturning at or near a critical level, and are points where strong turbulence is to be expected, as discussed by Ralph et 
al. [1997a], Thus, both in situ data and modeling results support the notion that this highly turbulent ER-2 flight was due to 
strong breaking of mountain waves at or near critical levels in the stratosphere. 

1.3 Fourier and Wavelet Spectra of DC-8 MMS Fluctuations 

Bacmeister et al. [1996, 1997] demonstrated that wavelet-based analysis of MMS data from accumulated ER-2 
flights was a useful and illuminating way of spectrally analyzing these data. Initial Fourier-based power spectral analysis of 
MMS data from SUCCESS was provided by Chan et al. [1998], and some preliminary wavelet-based spectral analysis was 
provided by Demoz et al. [1998], Both studies reported spectral shapes at small spatial scales ~k' 3 , where k is horizontal 
wavenumber, consistent with earlier stratospheric findings [Bacmeister et al., 1996], Evidence of a k' 5/3 turbulent regime was 
also noted at very short scales in some of the observations [e.g., Chan et al., 1998], Detailed wavelet and Fourier spectral 
analyses of data from all the flight segments have been undertaken. Results for one particular flight segment (segment 14) will 
be shown here, to give a feel for the sort of results such analysis has yielded. 

Figure 11 shows continuous wavelet spectra (right column) of unsmoothed (1Hz) zonal, meridional and vertical 
velocity fluctuations measured by MMS along flight segment 14 (see Figure 4), using the algorithm of Torrance and Campo 
[1998] with a Morlet wavelet basis (Gaussian-modulated sinusoids) and spectral oversampling. For zonal winds, we note that 
the wavelet spectrum identifies the 43 km mountain wave, but that it peaks most strongly in the eastward portion of the flight 
rather than in the westward part. This is consistent with the stronger isentropic oscillations and larger orographic falloffs in 
this region noted in Figure 3. Similar features are evident in the meridional velocity spectrum. Furthermore, the longer scale 
wave of wavelength -170 km also appears strongly in the wavelet spectra. Peaks are also evident at an intermediate scale -90 
km, which may be due to spectral leakage caused by a non-optimal choice of wavelet basis functions. 

Wavelet spectra of vertical velocities (bottom right panel of Figure 11) show evidence of these wavelengths. 
However, they also show evidence of intense short-scale perturbations. Indeed, we see clear evidence of strong oscillations of 
-15 km wavelength at the westward and eastward extremities of the flight segment that show out clearly in the wavelet 
spectrum. Closer inspection of the wavelet spectrum also reveals significant subsidiary peaks at wavelengths -1.5 km 
(isolated blue peaks in the spectrum at flight distances 3910km, 3990 km, etc.), which are not obvious from direct inspection 
of the time series. These features appear to be correlated to phases of the 15 km vertical velocity wave, and so could be 
coherent turbulent and/or nonlinear wavelike structures associated with instability/overtuming of the primary wave [e.g., 
Fritts et al., 1997], Another possibility is that the DC-8 may be sampling the wake vortex that it created earlier during the 
outbound flight segment #13 (an encounter with wake vortices on May 4, 1996 was reported by Chan et al. [1998]). Note that 


5 see the NRL/MWFM web site at http://uap-www.nrl.navy.mil/dynamics/html/mwforc.html 
see http://cloud 1 .arc.nasa.gov/success/daily_summary/Flight_reports/960502.er2.html 
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the locations of these peaks correlate well with peaks in TEDRs for flight segment 14 shown in Figures 4e and 5e. Wavelet 
analysis is ideal here, has it captures all of these strong nonstationary oscillations in various part of these time series. 

It is interesting to compare these results with the findings of conventional Fourier autospectra, shown in Figure 12. 
Note that the autospectrum of vertical velocity oscillations reveals peaks —15 km and —1.5 km, as in the wavelet spectrum, 
although the autospectrum goes not give the locations of this variability in the time series. The nonstationarity of the 
oscillations at these scales leads to noticeable spreading of the power spectral peaks in Figure 12. Light curves on the velocity 
autospectra in Figure 12 show k (solid curve) and k 3/3 shapes. We note for all spectra a noticeable transition, from k ' shapes 
at larger scales, consistent with gravity waves [Bacmeister et al., 1996], to k' 513 spectral shapes at short scales, consistent with 
inertial-range three-dimensional Kolmogorov turbulence [e.g., Chan et al., 1998]. Note in particular that the vertical velocity 
peak at -1.5 km wavelength occurs within this k~ 513 spectral regime, suggesting that these structures may be coherent turbulent 
structures of some kind generated by mountain wave dissipation (e.g., vortex rolls, Kelvin-Helmholtz instabilities, etc.). 

1.4 Generalized Parcel Advection Methodology for Modeling Gravity Wave Influences on 
Minor Constituents 


It is well known that gravity waves (including mountain waves) significantly influence cloud formation [e.g., Banta, 
1990] as well as contrail evolution and dispersion [e.g., Schumann et al., 1995; see also section 1.1.2], Yet the way in which 
gravity-wave dynamics interact with chemistry, microphysics, diffusion and larger-scale advection remains poorly understood 
at present [Schumann et al 1995]. 

As a formal approach to these issues, we have developed a parcel-based approach to the problem, based on the 
concept of quasi-adiabatic advection of air under the influence of gravity waves. Such an approach is often assumed to be 
valid for gravity waves, but has never been formally justified. Indeed, the approach differs significantly from traditional 
analyses of this problem for gravity waves, which involve linearized perturbation expansions of rate equations governing the 
chemical evolution of the relevant minor constituent. Such approaches require a background constituent profile that varies 
only gradually with height and horizontal distance, which is rarely the case for minor chemical species and aerosols in the 
upper troposphere [e.g., Newell et al., 1996]. Indeed, very thin aerosol and water/ice layers were also observed here [section 
1.1.2, Figures 6-8], The formal parcel-based analysis method we developed avoids such a restriction. Furthermore, the 
approach gels naturally with conventional parcel advection models of the larger-scale mean transport and quasi-horizontal 
diffusion of air by synoptic-scale motions [e.g., Brasseur et al., 1996], allowing interactions among these dynamics to be 
gauged more easily using gravity-wave-based additions to an existing modeling framework. 

Our method is based on first-principles theory of dry parcel advection. We derive results for a hierarchy of simplified 



Figure 13. Height-time surface renderings of vertical profiles of minor constituent number densities as a gravity wave passes through the 
layer, as simulated using the parcel advection methodology of Eckermann et al. [1998]. Plots show gravity wave interactions with: (a) a 
laminated ozone profile in the lower stratosphere; (b) a large narrow enhancement in mesospheric sodium densities. See Eckermann et al. 
[1998] in Appendix A for more details. 

problems which have been solved previously using conventional linearized perturbation expansions. The two methods give 
identical results. However, in environments where the minor constituent exhibits large spatial inhomogeneities (e.g., cloud 
and contrail distributions, tropospheric ozone etc.), the parcel method gives far superior results. This work was written up and 
accepted for publication in The Journal of the Atmospheric Sciences [Eckermann et al., 1998]. A reprint of the published 
article is provided in Appendix A. 

Sample results are shown in Figure 13. Panel a shows a gravity wave of vertical wavelength 2.3 km interacting with a 
laminated ozone profile in the lower stratosphere, which models the lidar observations of Gibson-Wilde et al. [1997], Similar 
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inhomogeneities (e.g., cloud and contrail distributions, tropospheric ozone etc.), the parcel method gives far superior 
results. This work was written up and accepted for publication in The Journal of the Atmospheric Sciences [Eckermann et 
al„ 1998]. A reprint of the published article is provided in Appendix A. 

Sample results are shown in Figure 13. Panel a shows a gravity wave of vertical wavelength 2.3 km interacting 
with a laminated ozone profile in the lower stratosphere, which models the lidar observations of Gibson-Wilde et al. [1997], 
Similar laminated ozone profiles occur in the troposphere [e.g., Newell et al., 1996; Langford et al., 1997]. We see in 
Figure 13a a periodic expansion and contraction of the width of the ozone lamina. Furthermore, Figure 13b shows that for 
narrower layers again, there is also significant vertical advection of the layer as well as contraction and dilation of its width. 
Standard linearized perturbation analyses cannot capture this behavior [Eckermann et al., 1998]. While the simulations in 
Figure 13 are nominally stratospheric and mesospheric, they assume tracer-like behavior (i.e. no radiation or chemistry) and 
so are directly applicable to the distortion and advection of analogous tropospheric layers, such as the narrow cloud and jet 
contrail layers sampled during SUCCESS. 

We have tested the method further in microphysical simulations of wave-induced effects on polar stratospheric 
clouds (PSCs). This is a sensible first extension of the method to non-tracers, since the chemistry and microphysics required 
to simulate PSCs is usually significantly simple- than for tropospheric clouds, thereby allowing more straightforward initial 
tests using the new modeling methodology. This was achieved by using the method to quantify wave-induced adiabatic 
cooling of air parcels by mountain wave upwelling, which was then used to extend NRL/MWFM to generate global maps of 
P°stcast mountain-wave-induced temperature perturbations [Eckermann and Bacmeister, 1998], The work has led to 
significant published work on PSC formation by stratospheric mountain waves in the Arctic and its possible influence on 
the ozone budget [Carslaw et al., 1998a, 1998b, 1999]. We also used the parcel method in a more theoretical assessment of 
how typical random spectra of gravity waves might influence ambient distributions of aerosol particles [Bacmeister et al., 
1998], All four papers have appeared in press or have been accepted for publication, but only the latter paper is included in 
Appendix B (to save space). These published works have come as spinoffs of the theoretical work initiated here as part of 
the current research project. Applications of the insights gained in these stratospheric simulati ons will be used to study the 
SUCCESS wave cloud and contrail data. Indeed, the findings in these stratospheric simulations may have quite direct 
relevance to upper tropospheric ice clouds and subvisible cirrus, since their microphysics and chemistry have man y 
similarities to PSCs [e.g., Jensen and Toon, 1994; Borrmann et al., 1996; Jensen et al., 1998]. The enhanced NRL/MWFM 
model will also play an important forecasting role in flight planning for future NASA aircraft campaigns, such as the SAGE 
III Ozone Loss and Validation Experiment (SOLVE), planned for the Arctic in the winter of 1999-2000. 

There was another unanticipated theoretical spinoff from this parcel-based work. In developing the parcel 
advection code, it became clear that superposing gravity waves would, in essence, “advect” parcels of constant potential 
temperature across other waves’ phase fronts. In other words, gravity waves should advect not only minor constituents, but 
the phase structures of other waves! Such a process (if valid) could lead to waveform distortions, which in turn should be 
relevant to measured profiles of wave-induced fluctuations and their spectra (e.g., section 1.3). Thus, a simple parcel-based 
advection model of superposing gravity waves was developed to assess these ideas. Waveform distortions did indeed arise. 
Furthermore, the spectra of the vertical profiles yielded regular quasi-universal AT 3 autospectra, where M is vertical 
wavenumber. Spectra such as these are routinely measured in the atmosphere and have been the subject of intense 
theoretical debate for more than a decade [e.g., Smith et al., 1987; Hines, 1991]. These findings were written up, and the 
paper has been accepted for publication in Geophysical Research Letters [Eckermann, 1998], A reprint of the published 
article is provided in Appendix C. These results provide new insights into mesoscale dynamical interactions and spectral 
formation. 

1.5 Numerical Model Simulations of the Mountain Wave 

In addition to NRL/MWFM simulations of mountain waves, we aim to 
simulate the waves directly using nonlinear numerical fluid dynamical models of 
flow over orography. The original goal was to work with Dr. J. T. Bacmeister and 
Dr. J. Rottman in using the MM5 model for this purpose. However, during this 
period both individuals have taken new jobs in other fields and are no longer be 
involved in the project. Presently, we are working with Dr. Andreas Dombrack in 
Germany on this problem, who uses the MM5 model to simulate mountain waves in 
the Arctic. This has proven useful, since Dr. Dombrack uses a modified version of 
MM5 that has been extended well into the stratosphere (lid at -27 km). This affords 
the potential to study not only the mountain wave activity intercepted by the DC-8 
near the tropopause, but also mountain wave-induced turbulence intercepted by the 
ER-2 in the stratosphere. 

To date, we have conducted some simplified studies using a two- 
dimensional nonhydrostatic numerical model of flow over orography (the 



Figure 14: Isentropes generated 
by eastward flow over Gaussian 
hill using the UNSW/NRL 
mountain wave model. 


Page 7 of 15 


NAS5-97247; Mesoscale Variability in SUCCESS Data (PI: S. D. Eckermann) 

Progress Report: Oct 1, 1997 - Sep 30, 1998 


2. Analysis 

2.1 Interpretation of Results Obtained to Date 


Dean-Day et al. [1998] presented an initial analysis of mountain wave activity inferred from the DC-8 MMS and 
MTP during this May 2 SUCCESS flight. We have expanded upon their preliminary analysis of these data. Combining the 
many types of data acquired aboard the DC-8 yields a great deal of information about the mountain waves. 

Flow on May 2, 1996 was from the west, corresponding to air flowing from left to right in Figure 3. From Figure 3, 
we also note that the underlying topography consists of a series of mountainous peaks, followed by a sharp drop in mean 

topographic elevation from 106°W to 105°W 
(the Continental Divide). Above this drop-off 
zone, strong isentropic dips are observed by 
the MTP in the upper troposphere, followed by 
sharp rises again further downstream. Very 
similar isentropic dips and rises have been 
observed in this region during the development 
of mountain-wave-induced downslope wind 
storms near Boulder. For example, Figure 15 
shows isentropes inferred from in situ aircraft 
measurements in this region during the January 
11, 1972 downslope windstorm event [Lilly 
and Zipser, 1972; Lilly, 1978], While the 
amplitudes in Figure 15 are much larger, the 
upper-level isentropic dip and its relation to 
the falloff zone of the Continental Divide are 
very similar to those seen in the MTP 
isentropes in Figure 3. Similar features are also 
observed in model simulations [e.g., Durran, 
1986]. Turbulence was also encountered to the 
right of the isentropic dips in Figure 15, due to 
overturning. Enhanced TEDRs were also found 

The larger-scale waves observed here remained quasi-stationary over a few hours, as opposed to the highly 
nonstationary mountain waves recently measured near Boulder and analyzed by Ralph et al. [1997a, 1997b], However, 
shorter scale transient waves were also observed in the MMS vertical velocities (see Figures 5c), with associated small-scale 
turbulence. Turbulent vertical velocities occasionally showed coherent structures, as evidenced by wavelet and Fourier 
spectral analysis (Figures 11 and 12), due either to breaking structures in the mountain wave or remnant wake vortices 
resampled by the DC-8. Spectral shapes were generally consistent with earlier stratospheric studies of airborne MMS data 
[Bacmeister et al., 1996], sometimes showing a clear transition from gravity wave to turbulence dynamics (Figure 12). 

Very thin “laminated” cloud aerosol and water/ice layers were observed just below the tropopause, supporting earlier 
observations of sharp upper tropospheric layering from the DC-8 [e.g., Newell et al., 1996] and from ground-based data [e.g., 
Langford et al., 1997], Upwelling by the mountain wave is identified in MTP data downstream of the mountain near flight leg 
D, where wave-cloud aerosol and water/ice are measured (Figures 6 and 7). This supports earlier microphysical simulations 
which suggested that mountain-wave-induced upwelling yielded sudden adiabatic cooling, which in turn led to heterogeneous 
freezing and wave cloud formation along flight leg D [Jensen et al., 1998].. 

2.2 Recommended Further Action 

The results discussed in this result raise a number of additional questions and avenues for further research. Some of 
the more immediate are: 

• further analysis of the May 2, 1996 data, specifically meteorological variability along other flight legs and correlations 
between measured constituents and winds/temperatures along the various flight tracks. 

• computation of more derived data analysis products (e.g., momentum fluxes, directionality, etc.) 

• similar analysis of data during the April 30, 1996 mountain wave and wave-cloud encounter over New Mexico. 

• further modeling of mountain wave activity observed in these data using existing model suites. 


s&mttmn ~ mim zm $ no s t mumm* mz 



Figure 15: Isentropes measured by aircraft on January 11, 1972, during a 
major downslope wind storm near Boulder, Colorado. Dotted lines show 
flight trajectories, plus signs show regions where significant turbulence 
was recorded onboard the aircraft (after Lilly [1978]). 

in similar locations for the mountain wave observed here (see Figure 5e). 
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• combined wavelet and Fourier analysis of data from other flights, to see whether the structures observed here using this 
method (e.g.. Figure 12) can yield information on other dynamics (e.g., wake vortices, turbulent wave-breaking 
structures, cloudy and clear-air turbulence, short-scale gravity waves, etc.) 

• further analysis of role of mesoscale variability in dispersion and advection of cloud and contrail layers. 

2.3 Relation to Ultimate Objectives of the Research Contract 

Significant progress has been made in the original objectives of the research project. Specifically: 

• in-depth data analysis revealed a wealth of information about mountain waves, turbulence, vortices and so forth during 
the flight of May 2, 1996, when wave clouds formed and were intercepted by the DC-8. 

• a formalized parcel advection methodology has been developed for modeling the way in which gravity waves interact 
with minor chemical constituents and narrow cloud/contrail layers (see Figure 13). This methodology provides a crucial 
theoretical framework within which more detailed studies can be formulated. Examples include condensation via wave- 
induced upwelling [e.g., Carslaw et al., 1998b; Jensen et al., 1998], mesoscale dispersion of contrails, interaction with 
large-scale dynamics and turbulence, and many other important dynamical unknowns [e.g., Schumann et al., 1995], 

• enhancements to the NRL/MWFM model to forecast/postcast temperature drops due to mountain wave upwelling 
[Eckermann and Bacmeister, 1998], which can now be used as inputs into global microphysical models of cloud and 
contrail formation [e.g., Carslaw et al., 1998a, 1998b, 1999], 

• progress in both Fourier-based and wavelet-based spectral analysis of airborne data has been made. Results in section 1.3 
indicate that careful analysis of results from both techniques can provide very useful dynamical information on waves, 
turbulence, wake vortices, wave breaking, and so forth. 

• comparison with model results has already occurred. NRL/MWFM simulations succeed in simulating the large mountain 
wave amplitudes measured on the DC-8, as well as the strong mountain wave-induced turbulence that buffeted the ER-2 

during the May 2, 1996 flights. Initial UNSW/NRL mountain wave simulations generate isentropes similar to 
observations. 

• science emanating from this research has already resulted in a significant body of publications in peer-reviewed journals 
[Carslaw et al., 1998a, 1998b, 1999; Bacmeister et al., 1998; Eckermann, 1998; Eckermann et al., 1998] 
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ABSTRACT 

Existing analytical models of wave-induced minor constituent fluctuations result from linearized perturbation 
expansions of a rate equation governing either number density or mixing ratio, whereas many numerical models 
now use isentropic parcel advection methods to simulate these effects. Exact relationships between the two 
approaches are not currently clear for gravity waves. Here, the parcel advection method is formalized and applied 
to derive analytical formulas for the response of vertical constituent profiles of arbitrary shape to adiabatic 
gravity wave displacements. These relations are compared to corresponding formulas from standard linearized 
perturbation analyses. Both methods accurately model perturbations produced by nondissipating hydrostatic 
gravity waves within idealized vertical tracer profiles. Both methods can also model wave-induced perturbations 
of minor constituents with shorter chemical lifetimes. This is demonstrated by using the parcel method to 
reproduce previous results for wave-induced fluctuations in upper- stratospheric ozone. The parcel-based approach 
yields more accurate models of nondissipating hydrostatic gravity wave effects on tracer distibutions with sharp 
spatial gradients. The authors demonstrate this by using the method to model wave-induced modulations of 
sporadic sodium layers in the mesosphere and of ozone laminae in the lower stratosphere. Conversely, the parcel 
method does not accurately model tracer density perturbations produced by nonhydrostatic waves, or situations 
in which the photochemical response of the constituent leads to significant diabatic feedback on the perturbing 
wave. & 


1. Introduction 

High-resolution observations of minor (trace) con- 
stituents in the middle and upper atmosphere often re- 
veal small-scale perturbations superimposed upon the 
mean background distribution. Such features can be pro- 
duced by any number and combination of dynamical, 
radiative, and chemical processes. For trace constituents 
with long lifetimes, however, dynamics must be pre- 
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dominately responsible for any observed variability of 
that constituent’s local number density over shorter 
timescales. 

One important source of such variability is waves. 
For instance, vertical soundings at middle and high lat- 
itudes often reveal large enhancements or depletions of 
lower-stratospheric ozone (a tracer at these heights) con- 
fined to narrow ( — 1—2 km) vertical layers (e.g., Reid 
and Vaughan 1991; Bird et al. 1997). These structures, 
known assortedly as filaments or laminae, result from 
complex synoptic-scale flow patterns produced by 
breaking planetary Rossby waves, which yield air par- 
cels at different stratospheric levels with vastly different 
histories and, hence, constituent densities (e.g., Newman 
and Schoeberl 1995; Orsolini et al. 1997). 

Oblique oscillatory displacements of air parcels from 
their equilibrium positions produce accompanying den- 
sity and temperature changes that can also perturb con- 
stituent profiles. Oscillations due to various types of 
planetary waves have been found to induce related per- 
turbations in ozone (e.g., Newman and Randel 1988; 
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Prata 1990; Hess 1990; Randel and Gille 1991; Randel 
1993; Stanford and Ziemke 1993; Wirth 1993; Ziemke 
and Stanford 1994; Engelen 1996) and other trace con- 
stituents (e.g., Randel 1990; Salby et al. 1990; Lim- 
pasuvan and Leovy 1995). On smaller scales, internal 
gravity waves can also give rise to fluctuations in con- 
stituent profiles (e.g., Reber et al. 1975; Gardner and 
Shelton 1985; Hedin and Mayr 1987; Sugiyama 1988; 
Hoegy et al. 1990; Wilson et al. 1991; Alexander and 
Pfister 1995; Langford et al. 1996; Bian et al. 1996; 
Bacmeister et al. 1997). 

In fact, since observational data on middle and upper 
atmosphere constituents are often better and more abun- 
dant than wind and temperature data, analyses of fluc- 
tuations in long-term measurements of constituent den- 
sities have been used to infer climatologies of both plan- 
etary waves and gravity waves (e.g., Randel and Gille 
1991; Senft and Gardner 1991; Randel 1993; Ziemke 
and Stanford 1994; Collins et al. 1996). However, since 
the response of constituents to waves is usually not 
straightforward, wave properties must be inferred from 
such observations by using a wave— tracer interaction 
model to convert the measured constituent response into 
a wave-related quantity. Such conversion methods have 
been derived to date by applying standard linearized 
perturbation methods of gravity wave theory to rate 
equations governing constituent densities or mixing ra- 
tios. Such derivations yield analytical formulas for 
translating conventional gravity wave oscillations into 
a corresponding response in the constituent profile (e.g., 
Dudis and Reber 1976; Chiu and Ching 1978; Gardner 
and Shelton 1985). Such relations can then be used to 
convert observed constituent fluctuations into wave-re- 
lated quantities (e.g., Gardner and Voelz 1985; Randel 
1990; Senft and Gardner 1991). 

Recently, detailed numerical models have provided 
new insights into the way minor constituents respond 
to planetary wave and gravity wave advection (e.g., 
Fritts et al. 1993, 1997; Waugh et al. 1994; Roble and 
Shepherd 1997). Rich complex structures in minor spe- 
cies distributions are simulated, which often resemble 
observations (e.g., Fritts et al. 1993, 1997) but are not 
predicted by standard formulas that convert between 
observed tracer variability and a plane monochromatic 
wave. Due to their complexity, however, these models 
cannot as yet be used for operational extraction of wave 
properties from observed constituent fluctuations in the 
atmosphere. On the other hand, simplified versions of 
these models, which retained the minimum terms needed 
to simulate the basic effects, might lead to improved 
analytical models of wave-induced perturbations of mi- 
nor constituent profiles. 

We investigate this possibility here by considering 
one class of these numerical models, the so-called tra- 
jectory or parcel advection models. They have proved 
successful in simulating finescale tracer structures in the 
stratosphere, often using wind and temperature data with 
much coarser spatiotemporal resolution (e.g., Newman 


and Schoeberl 1995; Orsolini et al. 1997). The models 
work well because their assumptions of adiabatic motion 
and long chemical lifetimes are well satisified for many 
constituents in the lower stratosphere, where diabatic 
heating and turbulent diffusivities are both small. How- 
ever, constituents with shorter chemical lifetimes can 
also be simulated by these models if a suitable parcel- 
based chemistry scheme is used (e.g., Lutman et al. 
1997). Even in the troposphere, where diffusion and 
diabatic heating can be large, parcel advection models 
that parameterize these effects have also proved to be 
useful tools (e.g., Kao et al. 1995; Rind and Lerner 
1996). 

Given that parcel advection models often perform 
well even after many days of integration, the same mod- 
els should also perform well if the synoptic-scale ad- 
vection patterns were supplemented with wave-induced 
advection effects since their timescales (periods) are of 
the order of days for planetary waves and less than a 
day for gravity waves. Furthermore, since many gravity 
waves and planetary waves propagate conservatively 
through the stratosphere, an assumption of adiabatic 
wave-induced advection should also be well satisfied. 
Indeed, limited tests of gravity wave effects in such 
models have used this sort of approach (Pierce et al. 
1994; Jensen and Thomas 1994; Carslaw et al. 1998). 

Thus, we formalize the basic physical concepts un- 
derpinning parcel advection models in section 2, assess 
the utility of the method in gravity wave studies, then 
apply it in section 3 to derive analytical responses of 
tracers to wave-induced displacements. We compare our 
results throughout with corresponding results using the 
traditional approach of expanding state parameters into 
the sum of a mean and perturbation term, then isolating 
and linearizing the perturbation terms in the governing 
equations. Results of these analyses and intercompari- 
sons are summarized in section 4. 

2. Air parcel approach 

This approach is based on the classic meteorological 
problem of vertically displacing a small dry parcel of 
air from its equilibrium height, where it is “well mixed” 
and in thermal equilibrium with the surrounding at- 
mosphere. We assume that the displacement is small 
and that the parcel responds adiabatically while always 
maintaining a pressure equal to that of the surrounding 
atmosphere, following standard conventions (see, e.g., 
section 2.5.1 of Wallace and Hobbs 1977; section 9.1 
of Iribarne and Godson 1981). 

a. Conserved quantities in the parcel 

1) Tracer mixing ratio 

We consider a constituent of number density n within 
an air parcel of total number density n M . For a trace 
constituent n/n M <C 1, any changes in n produce neg- 
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ligible changes in n M , whereupon the following conti- 
nuity equations apply: 

dn K s 

• U = 0 (la) 

dn _ 

— + nV • U = R, (lb) 

where U is the local velocity vector, d/dt = d/dt + U • V 
is a time derivative in the Lagrangian (parcel following) 
frame ( d/dt is the ground-based Eulerian time deriva- 
tive), and R is the net production/loss rate for the trace 
constituent’s number density n. Eliminating V • U be- 
tween (la) and (lb) yields (Lindzen and Goody 1965) 


d_( n \ _ dq _ R 
dt\n M ) dt n M ' 


( 2 ) 


where q is the mixing ratio of the trace constituent. 

For passive trace constituents, R ~ 0, whereupon (2) 
yields 



(3a) 


or 


n 

q — — = const (3b) 

within the parcel. Thus, we get the well-known result 
that passive trace constituents are also tracers since the 
movement of a region of unique constituent mixing ratio 
traces the Lagrangian motion of air parcels. 


2) Potential temperature 


Air parcels advected by an inviscid adiabatic flow 
also conserve potential temperature 0 so that 



(4) 


For a horizontally homogeneous atmosphere, the back- 
ground potential temperature 0 is given by Poisson’s 
relation 


*(z) = T(z)(^y ^ = T(z,z 0 ), (5) 

where z is height, y is the ratio of specific heats, T(z) 
is the background temperature, and p{z) is the_back- 
ground pressure (e.g., Turner 1973). From (5), ®(z) is 
also equal to the temperature T(z, z 0 ) that an air parcel 
attains on being transported adiabatically a distance Sz 
from its rest height z to a new height z 0 = z + &. 

3) Potential density 

We can also work in terms of potential density, D(z), 
which is defined for a horizontally homogeneous back- 
ground atmosphere as 


D(z) = nz)(^J y = fKz,Zo) (6) 

(e.g. Turner 1973), where p(z) = Mn M (z ) is the back- 
ground total air density at height z, n M (z ) is the back- 
ground total number density, and M is the mean mass 
of an air molecule, assumed constant (no diffusive sep- 
aration). From (6), D(z ) also equals the density p(z, z 0 ) 
that an air parcel attains on being transported adiabat- 
ically from a pressure, p(z) and density p(z), at its rest 
height z to a reference pressure p(z a ) at the displaced 
height z 0 = z + Sz. It then follows (e.g., Wallace and 
Hobbs 1977) that 


for parcels on an inviscid adiabatic flow, and thus the 
potential density of any parcel within such a flow is also 
conserved. 


4) Potential tracer density 

Analogously, we can also define a potential density, 
A, for a given tracer. If the constituent in question be- 
haves as a tracer at all heights and has negligible hor- 
izontal gradients in its distribution, then, as in (6), we 
can define a background potential tracer density profile 
given by 

A(z) = M lc V(z) = M, c n(z)(^^j * = M K n(z , z„), 

ViW/ ( 8 ) 

where M, c is the molecular mass of the trace constituent 
and v{z) is the background potential number density of 
the tracer. From (8), we see that A(z) also equals the 
mass density of tracer M lc n(z, z 0 ) that an air parcel con- 
tains after being transported adiabatically from a pres- 
sure p(z) and tracer density M, c n(z) at its rest height z 
to a reference pressure p(z 0 ) at the new height z 0 = z 
+ Sz. It then follows that 

cl A dv 

-J=* = 0 (9 > 

for parcels on an inviscid adiabatic flow. Since number 
densities are more commonly considered in tracer stud- 
ies, we shall work in terms of v rather than A in sub- 
sequent analysis. 

b. Perturbations due to adiabatic vertical 

displacements 

Next we consider local responses of each conserved 
term to a small vertical adiabatic displacement £ of an 
air parcel from its equilibrium position z x to a new height 

z 2 = z 1 + £ (10) 

given the standard background pressure profile 
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(a) Potential Temperature 


(c) Potential Density 





rv FlG j Schematlc diagram of response (light dotted curve) of (a) potential temperature ©, (b) absolute temperature T, (c) potential density 
D, and (d) absolute density p in an air parcel (shown as a cube) as it is adiabatically advected a small vertical distance £ from its equilibrium 
rest position z„ where it is in thermal equilibrium with the background atmosphere. Background profiles of each quantity are plotted with 
solid curves. J v 


p(z 2 ) = p(Zi) exp 

(n> 

where H(z, ) = RT(Zi)/g is the pressure scale height at 
height z 1; and the latter expression in (11) is sufficiently 
accurate_ for small 4 given a background temperature 
profile T(z ) that varies slowly with height. 

The formal parcel-based analysis of this problem is 
given in appendix A, where departures (perturbations) 
of parcel-based quantities from surrounding background 
values are derived. The basic procedure for computing 
these perturbations is depicted schematically in Fig. 1 
for both temperature and density. In each panel in Fig. 
1 , an air parcel has been advected from its equilibrium 
height z, to a new height, z 2 , and in each case the adi- 
abatic response is given by moving the parcel along the 
light dotted “response curve.” The top panels of Fig. 
1 show constructions for the potential temperature and 



potential density: since both remain constant within the 
air parcel, their response curves are vertical. Corre- 
sponding constructions in terms of the absolute tem- 
perature and density are shown in the lower panels. 
Their response curves are given by Tfz,, z 2 ) and p(z,, 
z 2 ) in (5)-(6): hereafter, it proves expedient to use the 
abbreviated notation T(z 2 ) and p(z 2 ) for these expres- 
sions since z, stays constant and the functional depen- 
dence on z, is implicit from (10). Background profiles 
are shown with thick solid curves on each plot. Per- 
turbations are given by the difference between the adi- 
abatically varying value within the parcel and the local 
background value at height z 2 , as shown in Fig. 1. Using 
the abbreviated notation introduced above, this yields 
perturbation terms @'(z 2 ), D’(z 2 ), T'(z 2 ), and p'(z 2 ). 

Simplified equations for these perturbation terms fol- 
low when the displacement £ is sufficiently small that 
variations in both the background profiles and the ad- 
iabatic response curves between z, and z 2 are approx- 
imately linear in Fig. 1. This holds so long as any cur- 
vature of these profiles occurs over vertical scales L » 
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£• Since the formal parcel relations in appendix A are 
height-varying exponentials, L can be conveniently 
characterized by the scale height H x of the exponential 
relation for quantity X, so that \£/H x \ C 1. When this 
condition is satisfied, the exponential parcel relations 
can be accurately approximated by retaining only the 
first two terms from their MacLaurin series expansions. 
For temperatures and densities, these truncated Mac- 
Laurin series (TMS) relations yield the instantaneous 
simplified perturbation relations 


P'fe) _ gfe) _ = (r(z 2 )\ 

p(z 2 ) d(z 2 ) \©( z 2 )/ 


_ t ^( Zl ) 

g 

where 


( 12 ) 


H D ( Zl ) = H( Zl ) 


y - 1 , dH( Zl ) 
y d z 


N 2 ( Zl ) 


(13) 


is the potential density scale height (see appendix A). 
Equation (12) holds for l£/// D l C 1. 


c. Range of applicability to gravity waves 

Next we consider a sinusoidal form for the vertical 
parcel displacement from z x : 

£ = £'(*, y> z l9 t) 


— I sin (kx 4- ly + m Zl — cot + <p), (14) 

where l is the peak vertical displacement amplitude, <p 
is a (constant) phase offset, ( k , /, m) is the wavenumber 
vector, and a> is the oscillation frequency of the parcel 
about its equilibrium height z x . Considering potential 
temperature only for the moment, (14) and (12) yield 


©'&_?> z 2 , 0 

©fe) 


~N\ Zl ) 


£'(*> y> z x , t). (15) 


Since_£ is small (l£/// D l < 1), replacing the variable 
term ©(z 2 ) with the constant value ®( Zl ) introduces 
negligible errors in (15). Then, to a good approximation, 
we can remove the explicit transport of parcels to z 2 so 
that we can reexpress the potential temperature pertur- 
bation (15) as a sinusoidal expression of the form (14) 
at a fixed height z x : 

©'(*, y, z x ,t) = ©(Zj) sin(kx + ly + m Zl - cot + <p\ (16) 


where 


©fei) 


~N 2 ( Zl ) 

g 


®(Zi)^ 


d Z 


(17) 



The relations (14) and (18) have the same form as 
the simplified polarization relations that govern a hy- 
drostatic gravity wave of intrinsic frequency co . This 
occurs because, in the linear hydrostatic limit, a gravity 
wave advects air parcels adiabatically (e.g., Fritts and 
Rastogi 1985), has a small displacement amplitude 
£ (I &H d \ C 1), and produces a negligible pressure per- 
turbation within displaced parcels (e.g., Gossard and 
Hooke 1975). 

To assess the degree and range of this correspondence 
more quantitatively, a standard non-Boussinesq linear- 
ized perturbation derivation of the acoustic-gravity 
wave equations is outlined in appendix B. There it is 
shown that the parcel-derived polarization relations (18) 
are an acceptable approximation to the full acoustic- 
gravity wave polarization relations when the normalized 
quantities \a(m , (o) I and \b(m, to)\ are <Cl [see Eqs. (BIO) 
and (B 11)]. When this limit is not satisfied, wave-in- 
duced pressure perturbations become nonnegligible [see 
Eq. (B13)], so the parcel approach (which assumes zero 
pressure perturbation) becomes inaccurate. Sample cal- 
culations in appendix B show that the differences be- 
tween the full acoustic-gravity wave polarization re- 
lations and the parcel-based approximations (18) are 
small for \ z = lirllml ^ 20 km (see Figs. Bib and Blc). 
Thus, we use this as an approximate upper-wavelength 
bound on the validity of the hydrostatic gravity wave 
approximation and of subsequent parcel-based deriva- 
tions of gravity wave-induced minor constituent vari- 
ability. 

3. Gravity wave perturbations of minor 

constituents 

We now apply the parcel-based approach to model 
gravity-wave-induced perturbations of minor constitu- 
ent profiles. We begin by considering some simple ex- 
amples and compare the parcel-based derivations (ap- 
pendix A) with corresponding equations derived by a 
conventional linearized perturbation analysis (appendix 
B). 


a . Tracer with linear q( Z ) profile 

The parcel-based approach to this problem is depicted 
schematically in Fig. 2a and yields the perturbation for- 
mula (A 19). Since the background profile is linear, then 
q(z 2 ) = q(z x ) + (&q/d Z )C , whereupon (A19) and the 
wave solution (14) yield 


Transport of parcels to z 2 is now implicit within (16) 
and (17). Similar relations to (17) follow from (12) and 
(14) for relative perturbations of potential density, den- 
sity, and temperature, with peak amplitudes interrelated 
as 


q\x y y, z 2 , t) = - y, z 1? t). (19) 

As for (15), if l is small, then the explicit transport of 
parcels to z 2 can be removed in order that (19) can be 
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(a) Tracer Mixing Ratio 



(b) Tracer Number Density 



(c) Potential Tracer Number Density 



Fig. 2. As in Fig. 1 but for (a) tracer mixing ratio q, (b) tracer 
number density n, and (c) potential tracer number density v of the 
parcel. 


reexpressed as a wave solution of the form (16), with 
peak amplitudes related as 


( 20 ) 

A standard linearized perturbation analysis of (3a) 
using the methodology in appendix B yields the lin- 
earized perturbation mixing ratio equation (e.g., Randel 
1990) 


— Lcoq + 



( 21 ) 


which, since w = — also yields (20) at a given 
height Z\. Thus the two approaches give identical results 
in this case. 


b. Tracer with a linear n(z) profile 

The parcel-based approach to this problem is depicted 
schematically in terms of tracer number density in Fig. 
2b and in terms of potential tracer number density in 
Fig. 2c. The mathematical derivations are set out in 
appendix A. Because n(z ) is linear and v(z) is approx- 
imately linear over the interval £ then TMS expansions 
of the parcel relations for tracer perturbations [e.g., 
(A25b)] are accurate, leading to the simplified TMS 
result 

F'(x, y, Za, t) = n'(x, y, z 2 , t) ^ £'(x, y, Zu t) 
v(z 2 ) n(z 2 ) H v (Zx) ( 

[see (A26)], where 




( 1 + 1 dn(z,) \ 

\yH(z t ) n(z,) dz ) 


(23) 


is the “potential tracer density scale height” by analogy 
to the potential total density scale height (13) [see (A7)]. 
Since n(z ) is linear, the accuracy of the TMS relation 
(22) is controlled only by the curvature of the expo- 
nential response curve (A21) shown in Fig. 2b and so 
is accurate for \C!yH\ <C 1. 

As for (15) and (19), if l is small, then the explicit 
transport of parcels to z 2 can be removed from (22) so 
that this equation can be reexpressed as a wave solution 
of the form (16) at the equilibrium height z 1; with peak 
amplitudes related as 


n(z i) 
n(Zi) 


1 

ytffe) 


+ 


1 dn 
n(z,) dz 


C, 


(24) 


where dn/dz is a constant here. 

Equation (24) also follows from a linearized pertur- 
bation analysis using the nonhydrostatic isothermal 
acoustic-gravity wave equations of Hines (1960) (Dudis 
and Reber 1976; Chiu and Ching 1978; Weinstock 1978; 
Gardner and Shelton 1985; Sugiyama 1988). A similar 
analysis in appendix B shows that the full expression 
for nonhydrostatic acoustic-gravity waves is [see (B20)] 


w(z i) _ _ 1 1 dn 

"&) \yH(zi) n(z,) dz 


jV 2 fa) 

8 


[a(m, (o) + ib(m, &>)] 


l 


(25) 


Therefore (24) is valid only when I a(m, <w)l and I b(m, 
«)l are negligible, which, as discussed in section lc and 
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Nonlinear n(z) Profile 



Fig. 3. As in Fig. 2b but illustrating the complications introduced 
by advecting a parcel across a nonlinear n(z ) profile. 


appendix B, holds only for hydrostatic gravity waves 
of \ z ^ 20 km. Thus, (24) is a hydrostatic result (see 
also Dudis and Reber 1976), despite the fact that non- 
hydrostatic wave equations are needed to arrive at this 
result in some linearized perturbation analyses (e.g., 
Gardner and Shelton 1985). Note that no such caveats 
applied to the mixing ratio relation (20) — it holds for 
both hydrostatic and nonhydrostatic waves. Therefore, 
mixing ratio is the more straightforward quantity to con- 
sider when converting between minor constituent per- 
turbations and displacement perturbations due to non- 
hydrostatic gravity waves or “tall” (fast) Kelvin waves. 

c. Tracer with a weakly nonlinear n(z) profile 

Of course tracer density profiles cannot be infinitely 
linear, but instead usually approximate some type of 
Chapman curve. If the curvature of the profile is gradual 
with height, then, for hydrostatic waves, (24) is often 
used to model wave-induced perturbations, with dn(zi)l 
dz reevaluated at each new height z x in (24), as in (23) 
(e.g., Chiu and Ching 1978; Gardner and Voelz 1985; 
Senft and Gardner 1991). For this to be accurate, the 
profile needs to be approximately linear only over the 
limited range z x ± as shown in Fig. 2b. However, 
this approximation may not always be sufficiently ac- 
curate for certain profile shapes, as shown schematically 
in Fig 3. Thus, tests of its adequacy for a given profile 
are required. 


In the parcel-based approach, the assessment is 
straightforward. For an exponential n{z) profile, for in- 
stance, (22) and hence (24) are accurate for \£/H v \ <C 
1, from (A6a) and (A6b). For arbitrary profile shapes, 
we simply avoid TMS approximations and use the gen- 
eral parcel solution (A21), which, using (14), yields the 
normalized perturbation relation 


n\x, y, z 2 , t) 

n(z 2 ) 


n (Zi) (-Z'(x, y, z„ t) 

w exp r~^r 


1 . ( 26 ) 


from (A23b), where z 2 = z, + £'(x, y, z„ t). Note that 
(26) yields a nonsinusoidal solution. In terms of poten- 
tial tracer density v, this is equivalent to using (A6a) 
or (A8) instead of the approximate relation (A6b). Direct 
comparison of results using (24) and (26) then indicate 
whether the former equation is sufficiently accurate for 
modeling wave perturbations of a given n(z) profile. 

Such assessments are less straightforward within the 
framework of a linearized perturbation analysis (e.g., 
Chiu and Ching 1978). Gardner and Shelton (1985) de- 
veloped an iterative method for deriving higher-order 
corrections to (24) due to the curvature of a given n(z) 
profile, using perturbation expansions based on an as- 
sumed form for the solution. Their most general solution 
took the form 


n(x, y, z, t) 



z - yH(z) In 


1 + 


A'(x, y, z, t) 

y ~ 1 


1 + A'(x, y, z , t) 

y ~ 1 


(27) 


where A'(x, y, z. t ) is a normalized gravity wave— induced 
perturbation, which, for small-amplitude waves in an 
isothermal atmosphere (Hines 1960), is approximately 
equal to the relative density perturbation of the wave 
[see Eq. (38) of Gardner and Shelton 1985], It then 
follows that (27) corresponds approximately to the par- 
cel solution (A21) derived in appendix A since 1 + 
A'(x, y, z, t)/(y — 1) = 1 + t' (.v, y, z. I)lyll(z). which 
in turn is a TMS expansion of exp [£'(x, y, z, t)lyH(z)]. 

Gardner and Shelton converted their solutions into 
perturbation formulas at a set height z, using MacLaurin 
series expansions. They derived the following extension 
of (24) in which the effects of both first- and second- 
order vertical derivatives of n(z) were retained: 


n'(x,y,z l ,t ) -1 

1 + 1 

( r^ r fei)\3«(z 1 ) 

p'(x, y, z u t) 

i 

1 

i, «(Zi) ) dz 

P(Z i) 


+ 


1 

(7 ~ D 2 


1 + 


{ 3yH(z,) \ dn (z,) + 
V 2n(z,) ) dz 


(y 2 H>(Zi)\ 

1 d 2 n(Zi) 

(p\x, y, Zi, t)\ 

\ 2 7f(Zi) ) 

dz 2 

\ P(z.) ) 


(28) 
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where p'(x y y, z l5 t) is the gravity wave density oscil- 
lation, which in this case is related to l through (14) 
and (18) using isothermal relations ( dH(z x )!dz = 0) to 
evaluate A 2 (z,) in (13). Like (26), Eq. (28) is also a 
nonsinusoidal solution. 

In Fig. 4, we compare (28) with the linear profile 
solution (24) and the parcel solution (26). Since sodium 
acts as a passive tracer to gravity waves above —85 km 
(e.g., Hickey and Plane 1995), we use an n(z) profile 
based on a Gaussian model of the mesospheric sodium 
layer, similar to the one considered by Gardner and Shel- 
ton (1985). The profile is plotted as the thick light curve 
in Figs. 4b and 4e. We perturb this profile with a hy- 
drostatic gravity wave (14) of vertical wavelength \ z = 
10 km, t = 800 m (both representative mesospheric 
values), and constant phase (p 0 . Sample waves are shown 
in Figs. 4a and 4d as a function of height z at a given 
location (x 0 , y 0 ) at two different times t. The corre- 
sponding perturbations to the background tracer profile 
are plotted in Figs. 4b and 4e, and the normalized per- 
turbations are isolated and plotted in Figs. 4c and 4f for 
the parcel solution (26) and the perturbation formulas 
(24) and (28). The parcel solution (26) is also replotted 
with a dotted curve in Figs. 4c and 4f at the undisturbed 
height z = Zi, instead of at the perturbed height z = z 2 
= Zi +'i'(x 09 y 0 , Zi, t) (solid curve), to help separate 
the deviations that arise due to different physical lo- 
cations for each parcel (i.e., differences in z 2 and z x ) 
and those due to incomplete characterization of the ef- 
fects of mean profile curvature. 

We see in Figs. 4c and 4f that the simple relation (24) 
is reasonably accurate here but that the higher-order 
correction (28) brings the solution closer to the parcel 
solution (26). The results in Fig. 4 also show many well- 
known features of wave-induced fluctuations of tracer 
number density profiles; for example, that the wave re- 
sponse is always larger on the underside of the layer 
than on the topside (e.g., Chiu and Ching 1978; Wein- 
stock 1978). This feature is easily appreciated in a parcel 
analysis by noting that all points on the unperturbed 
Gaussian profile in Figs. 4b and 4e follow the “adiabatic 
rarefaction” response curve (A21) when vertically ad- 
vected, as shown in Fig. 2b. This clearly yields greater 
deviations from the background profile on the bottom 
side of the layer than on the top side. Also evident in 
Figs. 4c and 4f is the well-known half- wavelength phase 
flip across the layer peak (Gardner and Shelton 1985). 

d. Tracer with a strongly nonlinear n(z) profile 

Tracer layers in the middle atmosphere often deviate 
significantly from the smoothly curved profiles consid- 


ered in Fig. 4. For example, the mesospheric sodium 
layer sometimes develops a large narrow spike of en- 
hanced Na density, known as a sudden (or sporadic) 
sodium layer (SSL) (Clemesha et al. 1980; Kwon et al. 
1988; Nagasawa and Abo 1995). Similarly, ozone pro- 
files in the lower stratosphere near the edge of the polar 
vortex often contain large narrow enhancements or de- 
pletions, known as “laminae” (e.g., Reid and Vaughan 
1991). Neither the linearized perturbation relation (24) 
nor its second-order correction (28) are sufficiently ac- 
curate for modeling wave-induced effects within such 
profiles. 

However, parcel-based methods can be used since 
they require no prior assumptions about the shape of 
the background vertical profile. To illustrate this, we 
simulate an SSL by adding a narrow Gaussian layer to 
the original Gaussian in Fig. 4, choosing its parameters 
to resemble the “typical” SSL considered by Cox et al. 
(1993). The resultant profile is the thick light curve 
shown in each panel of Fig. 5. We perturb this profile 
with the same wave as in Fig. 4. Solid curves in Fig. 
5 show the perturbed profile as calculated using the 
parcel solution (A21) and the wave oscillation (14), with 
each panel showing results at six successive equispaced 
time values t spanning a full wave period lirla). We see 
that the SSL is not only advected up and down by the 
wave, but is also periodically expanded (Fig. 5b) and 
compressed (Fig. 5e). 

Detailed features of the response are more easily seen 
in surface renderings of successive profiles, shown in 
Fig. 6. Figure 6a shows the periodic oscillation of the 
SSL peak and the downward movement of wave phase 
fronts above and below the SSL, the latter associated 
with an approximately linear wave response within the 
basic Gaussian profile (Figs. 4c and 4f). The SSL is 
focused upon in Fig. 6b, which shows a sinusoidal ad- 
vection of its peak value and periodic contraction and 
expansion of its width. 

Figure 7 shows a representative lower- stratospheric 
ozone profile with a lamina superimposed (thick light 
curve), based on the lidar observations of Gibson-Wilde 
et al. (1997). Like sodium, ozone has a long chemical 
lifetime at these heights and therefore acts as a tracer 
to gravity wave motions. We perturb the profile with a 
gravity wave of £ = 200 m and \ z = 2.3 km, choices 
based on the gravity wave detected and characterized 
by Gibson-Wilde et al. (1997) at -12-18 km during 
their ozone measurement. Parcel results using (A21) for 
a full wave cycle are shown with solid curves in Fig. 
7. Again, we see not only vertical advection of the 1am- 


Fig. 4. Perturbations of a Gaussian number density profile n{z) [thick light curve in (b) and (e)] are plotted in (b) and (e), as produced 
by the gravity wave shown above it in (a) and (d), respectively. Corresponding relative perturbations are plotted below in (c) and (f ). Results 
using parcel advection methods (solid and dotted curves) are compared to the results (24) (dashed curve) and (28) (broken curve) based on 
linearized perturbation j>olutions. The thick light curves in (c) and (f) show the asymmetric amplitude envelope of the parcel solution, given 
from (A23a) by [n(z)/n(z ± £)] exp( + £/y H) - 1. An isothermal atmosphere of H = 7 km was used. 
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(a) cot = 0° 



tracer concentration (m' 3 ) 



tracer concentration (m* 3 ) 


(b) cot = 60° 



tracer concentration (m' 3 ) 


(e) cot = 240° 



tracer concentration (m 3 ) 


(c) cot = 120° 



tracer concentration (m' 3 ) 


(f) cot = 300° 



tracer concentration (m* 3 ) 


Fig. 5. Sequence of six profiles (at equispaced times t spanning one wave period) of a wave-perturbed Na layer (solid curve) as in Fig 4 
but with an SSL inserted at 95 km (based on Fig. 1 of Cox et al. 1993). The thick light curve shows the unperturbed profile in each case.’ 


ina, but compression and expansion of its width as well, 
features that are more obvious in the corresponding sur- 
face plot in Fig. 8. Effects such as these may account 
for some of the standard deviation in mean magnitudes 
and widths of ozone laminae found in collated ozone- 
sonde data (e.g., Fig. 7 of Reid and Vaughan 1991), as 
well as variability of other narrow atmospheric layers: 
for example, tropospheric chemical layers (Newell et al. 
1996), cloud structures (e.g., Carslaw et al. 1998), polar 
mesosphere summer echoes (e.g., Fig. 1 of Rtister et al. 
1996), and stratospheric aerosols (Shibata et al. 1994; 
Hansen and Hoppe 1997). Figure 8 also shows the reg- 
ular downward movement of wave phase fronts above 
and below the lamina, where the response approximates 
a linear wave solution of the form (24). 

e. Three-dimensional gravity wave advection of a 

tracer 

The features that produce laminated ozone profiles 
tend to occur in long narrow streamers of meandering 
air, which lead to large background ozone gradients not 
only vertically, but also horizontally (e.g., Waugh et al. 


1994; Orsolini et al. 1997). These gradients can be im- 
portant when attempting to isolate wave-related struc- 
tures in ozone measured within such environments (e.g., 
Danielsen et al. 1991; Teitelbaum et al. 1996; Bac- 
meister et al. 1997; Gibson- Wilde et al. 1997). Again, 
the approximate relations (24) and (28) cannot accu- 
rately simulate the effects of gravity waves on such 
three-dimensional tracer distributions. 

The parcel-based theory in appendix A assumed no 
horizontal gradients in the background tracer profile and 
cannot be directly applied here. However, the analysis 
is easily extended to accommodate full three-dimen- 
sional wave-induced displacements [*'(*,, y i9 z i9 t ), 
w'(x i9 y ,•» Zj, 0, C(x n y i9 Zj , 0] of a collection of parcels, 
i (i = 1, 2, 3, • • • ) at equilibrium positions (x i9 y i9 z,) 
within an arbitrary three-dimensional tracer distribution 
y t z). 

The mixing ratio response (A 19) generalizes to 
4l x i + X\x i9 y i9 z„ t) 9 y ; - + 7]'(x i9 y i9 z i9 t) 9 
Z,. + C(x i9 y i9 z„ t)] = q(x i9 y i9 z,), (29) 

which, as discussed in section 3b, holds for hydrostatic 
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and nonhydrostatic waves. If we can assume that back- responding extension of (A21) yields the number den- 
ground isentropic surfaces are horizontal, then a cor- sity response 


«[*, + *'(*„ y„ z„ t), y, + r,'(x„ y„ z „ t), z, + C(x„ y„ z,, /)] = n(x„ y„ z,) exp 


~C(Xj, y-„ z„ t) 
yH(z,) 


(30) 


for a background number density distribution n(x, y, z). 

As discussed in section 3b, (30) is accurate for hy- 
drostatic gravity waves only. In environments where 
large narrow horizontal gradients in background poten- 
tial temperatures/densities also arise, additional com- 
pression and expansion occurs as parcels are horizon- 
tally advected by the wave. In this situation, not only 
must the exponential in (30) be modified to account for 
this extra horizontal compression/expansion, but ex- 
tended wave equations may be needed that incorporate 
concomitant modifications to the wave motion (e.g., 
Lamb and Shore 1992). Since background isentropes 
rarely tilt from the horizontal by more than a few de- 
grees, then (30) should adequately describe gravity 
wave-induced perturbations of a three-dimensional dis- 
tribution of tracer number densities in most circum- 
stances. 

It should be noted that the parcel relations (29) and 
(30) are Lagrangian expressions. A typical application, 
then, would be in supplementing synoptic-scale parcel 
advection models with oscillatory gravity wave advec- 
tion effects (e.g., Pierce et al. 1994). The Eulerian tracer 
perturbations that such simulations will produce (e.g., 
fluctuations in vertical tracer profiles) will be quite com- 
plicated in general. 


f Gravity wave perturbations of shorter-lived 
constituents 


The parcel-based approach can also be used for short- 
er-lived constituents. We illustrate this by using it to 
reproduce some well-known results for ozone. 

While the lower- stratospheric ozone in Figs. 7-8 acts 
as a tracer, the ozone production/loss term in (lb), 
Ro 3 , becomes nonnegligible above —25-30 km. Assum- 
ing Chapman photochemistry, then 


*03 2 Jo 2 n o 2 


k 02 n M n 0i ’ 


( 31 ) 


where n x is the number density of constituent x, J x is 
its photolysis rate, and k x is its recombination rate (Bras- 
seur and Solomon 1984). As shown in appendix C, ap- 
plying a parcel-based analysis to (2) and (31) yields a 
linearized response of the photochemical term R Q3 to an 
adiabatic vertical displacement of the parcel of the form 


= -B(z i )q' (H (z 2 ) ~ C (z 2 )T'(z 2 ), (32) 

for small vertical displacements £ where <7 0 ,(z 2 ) = 
<7 o 3 (Z2) + <7o ,(z 2 ) from (A19) and T’(z 2 ) is the TMS 
temperature perturbation (A 13). Since R Q3 (z) = 0, then, 
for a linear q(z) profile. 


dt 


dq 

dz 


I w(z 2 ), 


(33) 


where w(z 2 ) = dUdt is the vertical velocity of the parcel 
at z 2 (see Fig. 2a). Thus (32) and (33) yield 


dq’oMi) 

dt 


+ w(z 2 ) 



-B(z 2 )qo 3 (z 2 ) 


~ C(z 2 )T'(z 2 ). (34) 

This equation was first derived by Lindzen and Goody 
(1965) using a linearized perturbation expansion of (2) 
and (31). 

We now substitute the wave solution (14). If £is small 
and the vertical variation of background quantities in 
(34) is gradual (as assumed in appendix C) and ap- 
proximately linear over the interval z, ± £ then explicit 
parcel transport to z 2 can be removed so that the lin- 
earized mixing ratio and temperature perturbations in 
(34) can also be expressed as wave solutions of the form 
(14) and (16). Substitution of those solutions into (34) 
yields, with the help of (21), (A5), and (13), the peak 
amplitude relation 



<7o 3 (Zi) 


JdaM^/dz C(Z,) ' 

dT( Zl )/d Z + r a L oj 


x T( Zl ), (35) 

where T a is the dry-adiabatic lapse rate. This is the basic 
ozone-temperature relation derived for Kelvin waves 
by Randel (1990). It illustrates the_two limiting cases 
ofthe response. For u 2 » B 1 and (C/w) 2 <C [(dq 0) /dz)/ 
(■ dT/dz + r„)] 2 , terms on the right of (34) are small 
while terms on the left are large and dominate, so (35) 
approximates the tracer solution (20). For and 

(C/o)) 2 [(Bq 0 Jdz)/(dT/dz + I),)] 2 , photochemistry 
dominates so that terms on the right-hand side of (34) 
are much larger than the parcel transport terms on the 
left. In this latter case, (35) tends toward the photo- 
chemistry-dominant limit 
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Fig. 6. Time-height surface rendering of the perturbed Na number 
density response in Fig. 5. Both plots span one wave period. Panel 

(a) shows variations over the full height range of the Na layer, while 

(b) focuses on the SSL. 


< 36 > 

which also follows on setting the left-hand side of (34) 
to zero. 

This parcel-based analysis assumes adiabatic parcel 
motion, whereas the perturbation q^(z 2 ) also produces 
a diabatic heating/cooling perturbation in the potential 
temperature equation. However, Zhu and Holton (1986) 
show that this leads to very minor departures from adi- 
abaticity for hydrostatic gravity waves. Thus parcel- 
based analyses also work well in simulating (hydro- 
static) gravity wave— induced perturbations of constit- 
uents with short chemical lifetimes, so long as the per- 
turbed chemical response does not lead to significant 
diabatic damping of the wave; otherwise, a linearized 


perturbation analysis should be pursued (e.g., Zhu and 
Holton 1986). 


4. Summary and conclusions 

This study has formalized the parcel advection meth- 
od, applied it to gravity waves, and compared its find- 
ings with standard results from linearized perturbation 
expansions of the governing equations. The combined 
results of section 3 have shown that the two approaches 
are equivalent in most circumstances. However, we also 
identified certain situations in which one or the other 
method gave more accurate results. 

The parcel-based method becomes inaccurate when 
simulating perturbations of tracer density profiles by 
nonhydrostatic gravity waves. For shorter-lived con- 
stituents, the parcel method also becomes inaccurate if 
the wave-induced photochemical response produces a 
significant diabatic feedback on the wave. In these in- 
stances, the results from traditional linearized pertur- 
bation analysis should be used instead. 

On the other hand, the parcel-based approach is the 
better method for modeling hydrostatic gravity wave 
perturbations of tracers with sharp spatial gradients in 
their background distribution. We illustrated this by us- 
ing it to model wave-induced modulations of sporadic 
sodium layers and ozone laminae (Figs. 5-8). 

While we have restricted attention here to gravity 
waves, the parcel methodology described here can also 
be applied to other types of waves, such as free and 
forced Rossby waves, mixed Rossby-gravity waves, 
Kelvin waves, and tides. For some of these waves, both 
meridional and vertical wave displacements are neces- 
sary to simulate their effects on tracers (e.g., Randel 
1993; Stanford and Ziemke 1993). Parcel advection 
models have already been used extensively to simulate 
the effects of stationary Rossby waves on stratospheric 
tracer distributions. 
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(a) cot = 0° 



(b) cot = 60° 



(c) cot = 120° 



(d) cot = 180° 



(e) cot = 240° 



(f) cot = 300° 



Fig. 7. As in Fig. 5 but using a representation of the ozone layer based on the data of Gibson-Wilde et al. (1997), with a lamina at 14 
km. Thick light curves show the unperturbed profile. Solid curves show perturbed profiles due to a gravity wave of A = 2 3 km and l = 
200 m. z J b 


APPENDIX A 



Fig. 8. As in Fig. 6a but showing the ozone number density re- 
sponse from Fig. 7. 


Perturbation Formulas Using the Parcel 
Methodology 

a. Background profiles 

The ideal gas law and the altimeter equation (1 1) yield 
the standard background density relation 

?fe) ' exp (sfc)) ” ?<Z|) exp fefe))' 

(Al) 

where 

H p ( Zl ) = ff( Zl )^l + ' (A2) 

is the density scale height. On invoking (5) and (6) at 
2 ! and z 2 , then, with the aid of (11), (Al) and (A2), 
background potential temperature varies with height as 
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= Qfe,)ffi;exp( ^ 2 Z ‘ )(y 1} 

Ttei) \ yH(z x ) 


(A3a) 


~ ®(z,) exp 


H D tei)' 


(A3b) 


while background potential density varies with height 
as 


Z>(z 2 ) = Z>(z,)^ex p^ 2 2,) 


= D(z,) exp 


p(z>) \ 


H d(z i)/’ 


where 


// D (z,) = //(z,) 


y - 1 + 

y dz 


W 2 (z,) 


(A4a) 

(A4b) 

(A5) 


is the potential density scale height, g is gravitational 
acceleration, and N(z,) is the background Brunt-Vaisala 
frequency at height z, . 

From (8), analogous relations for the background po- 
tential number density of a given passive trace constit- 
uent are 


V(Z 2 ) = V(Zy) 


n(z 2 ) 

«(Zi) 


exp 


f (z 2 ~ z,) \ 
l y »(z l ) ) 


(A6a) 


where 


= v(Zi) exp 



(A6b) 


H v ( Zl ) = -H( Zl ) 



H(Zi) (In (z,) \ 1 

dz ) 


(A7) 


is the potential tracer density scale height. The simpli- 
fied relation (A6b) is accurate when height variations 
of the TT(z , ) profile are approximately exponential over 
the interval £ For small £ (A6b) is a good approxi- 
mation whenever 7f(z) varies over a typical vertical 
length scale L » £ However, while background tem- 
peratures in (Al) and background densities in (A4a) 
usually obey such conditions, n(z ) profiles in (A6a) of- 
ten do not (see sections 3c and 3d), whereupon (A6a) 
or the integrated expression 


v(z 2 ) = v(Zi) exp 



+f 


-dZ \ 

H P (Z)J 


(A8) 


must be used. Note that the reformulation of (A6a) into 
a fully exponential form (A6b) and (A8) is similar to 
the approach taken by Lindzen (1981) in parameterizing 
the variation with height of WKB gravity wave equa- 
tions [see his Eqs. (11) and (12)]. 


b. Temperature perturbations 

Figure la constructs the perturbation of potential tem- 
perature produced by the adiabatic vertical displacement 
£ of an air parcel from its equilibrium position z,, as 
discussed in section 2b. Since © in the parcel is con- 
served, the displacement produces a perturbation from 
the background profile at height z 2 of 

©'(z 2 ) = ©(z,) - W(Z 2 ). (A9) 

From (A3b) and (A5), this implies a normalized per- 
turbation: 


Q'fe) 

®(z 2 ) 



■ -c _ -jV 2 (z.) 

Hotel) g 


(A 10a) 


_1 a©(z,) » 

©(z,) dz 


(A 10b) 


where (A 10b) follows from retaining the first two terms 
of the MacLaurin series expansion of the exponential 
in (AlOa) and is accurate for small £ (i.e., I£/// D (z,)l <c 
1). This simplification is referred to hereafter as a trun- 
cated MacLaurin series (TMS) expansion and from 
(A 10b) is clearly accurate when the background profile 
©(z) is approximately linear between z, and z 2 . 

The absolute temperature response is depicted in Fig. 
lb. Using (5) and (11), the temperature of the displaced 
parcel becomes 


7-fe) = 7-fe) «p( : ^y 2 ) (Alla) 
^Ttet)-^-, (Al lb) 

c p 

where C p is the mass specific heat at constant pressure. 
Equation (Allb) follows from a TMS expansion of 
(Alla) and reproduces the familiar result that the parcel 
temperature varies according to the (dry) adiabatic lapse 
rate F ( = g!C p for small £ This yields a temperature 
perturbation 

T'tet) = T(z 2 ) ~ T(z 2 ) (A 12) 

and a relative perturbation, according to (A5), (Alla), 
and (A 12), of 


T(Z 2 ) 


Ttei) 

T(z 2 ) 


exp 


( ~&y ~ D \ 

\ y Htei) ) 


- l 



- l 


.. ~C -^ 2 (z.) 

Hotel) g 6 ’ 


(A13) 


using a TMS expansion of the exponential. Thus the 
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relative perturbations of absolute and potential temper- 
ature are equal for small £ 


c. Density perturbations 

Figure lc constructs the perturbation of potential den- 
sity produced by the same adiabatic vertical displace- 
ment £ of an air parcel from its equilibrium position Zl . 
Since D in the parcel is conserved, then, as in (A9), 

D'(z 2 ) = D( Zl ) - D{z 2 ). (A 14) 

From (A4b) and (A5), (A14) implies a normalized per- 
turbation, 


D\z 2 ) ( £ \ . 

= exp fei) " 1 

_ . £ A^fa) 

H d (z x ) g 


(A15a) 

(A15b) 


tf'fe) = q(z 2 ) ~ q(z 2 ) 
= q(Zi) ~ q(z 2 ) 
and 


g'fa) = q(z t ) _ 
q(z 2 ) q(z 2 ) 

Since p - Mn M , then (3b) and (A 16) imply 
n(z 2 ) = q{Zi)n M {z 2 ) 


(A19) 


(A20) 


= <A21) 

where n(z ) is the background number density profile of 
the tracer. Thus, like the total density p in (A 16), the 
tracer number density n responds to adiabatic vertical 
displacements according to its adiabatic rarefaction 
curve (A21). As shown in Fig. 2b, this produces a per- 
turbation of tracer number density 


where (A15b) follows from a TMS expansion of (A15a). 

Figure Id shows the corresponding construction for 
the density perturbation. Using (6) and (1 1), the density 
of the adiabatically displaced parcel responds according 
to 


n'(z 2 ) = n{z 2 ) - n(z 2 ), (A22) 

whose relative magnitude is given by 

= wfe) _ j 
n(z 2 ) n(z 2 ) 


P(z 2 ) = p(z,) exp 


~£ 

yH(z t ) 


(A 16) 


"(Zi) ( ~£ \ 

n{z 2 ) eXP \y//(z,)j 


- 1 


which, by analogy to the adiabatic lapse rate for tem- 
perature (A lib), we refer to here as the adiabatic rare- 
faction curve. This yields a density perturbation 


W(z 2 ) 



- 1 


(A23a) 

(A23b) 


p'(z 2 ) = p(z 2 ) ~ p(z 2 ), (A 17) 

which, with the aid of (Al), (A2), (A5), and (A 16), 
yields to a relative density perturbation of 

p\z 2 ) = p(z 2 ) _ 
p(z 2 ) p(z 2 ) 



,, i N 2 ( Zl ) 

H d (Z\) g 


using a TMS expansion of (A 18a). 


(A 18a) 
(A 18b) 


d. Tracer perturbations 

Figure 2a constructs the perturbation mixing ratio q' 
induced by the adiabatic vertical displacement £ of a 
parcel from its equilibrium height Zl in an atmosphere 
containing a tracer with a background vertical profile 
q( Z ). Since q in the parcel is conserved according to 
(3), then 


= g'fe) + p'fe) + / VfeA / p'feA 

q(z 2 ) p(z 2 ) \q(z 2 ))\p(z 2 ))’ 


(A23c) 


where (A23c) is derived from (A23b) using (A 18a) and 
(A20). 

An equivalent derivation in terms of potential tracer 
number density is shown in Fig. 2c. Since v in the parcel 
is conserved, then 


v> (z 2 ) = Kzi) - v(z 2 ), (A24) 

which, with the aid of (8), (11), and (A23a), yields 


^'fe) = v(z 
V{z 2 ) v(z 2 ) 


(A25a) 


n{z i) 
n(z 2 ) 


exp 


-£ 

yH(z t ) 


- l 


(A25b) 


« , fa) 

n(z 2 )‘ 


(A25c) 


If £ is small and v(z) variations are of much longer 
vertical scale than £ then (A25c) can be simplified using 
a TMS expansion of (A6b), along with (A7) and (A25a), 
to yield 
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(b) Errors in magnitude of parcel formulae 


6P 


17771 error: 

A 

Irl error: 


2P 


10 
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K (km) 


40 


50 


(c) Errors in phase of parcel formulae 


A — 

T/T phase error. 


r phase error. 


10 


20 30 

\ (km) 


40 


50 


Fig. Bl. Vertical wavelength variation of a max (m) and 6 max (ra) (a), 
and deviations in magnitude (b) and phase (c) between the acoustic- 
gravity wave polarization relation (appendix B) and the corresponding 
parcel-based approximation (appendix A) relating l to perturbations 
in relative temperature (dashed line) and relative density (dotted line). 


= n'(z 2 ) ^ C 

v{z 2 ) n(z 2 ) //„(z,)' 


(A26) 


The same expression also results from a direct TMS 
expansion of (A23a). Note the similar forms of (A26) 
and (A 15b). 


the following perturbation equations govern a zonally 
propagating nonhydrostatic gravity wave in a compress- 
ible rotating unsheared inviscid atmosphere (e.g., Gos- 
sard and Hooke 1975): 

-Lcou - fv = -i i:4>, (B3) 

-uov + fu = 0, (B4) 

-LMW + rg = —(iM - (B5) 

(B6) 

-L(or + iku + ( iM - H~')w = 0, (B7) 

where c s is background sound speed, ( u , v, w) is the 
peak velocity amplitude, r = p/p is the relative density 
amplitude, and d* = p/p, in which p is the pressure 
amplitude. The full nonhydrostatic dispersion relation 
for gravity waves [e.g., (23-7) of Gossard and Hooke 
(1975)] follows from the equation set (B3)-(B7); that is, 


m 2 = k 2 


/ V 2 - aA 

U 2 - f 2 ) 



(B8) 


Polarization relations follow from (B3)-(B7). Elimi- 
nating 4> between (B5) and (B6) and noting that w = 
— toil, then 


r = 


£ 

P 



[1 + a{m, (o) + i b(m, &>)]^, 


(B9) 


a(m, (o) = (1 - (oVN 2 ) 


(g\ 

-r 

\c 2 
\ s J 

m 2 + r 2 


(B10) 


b(m, (o) = (1 - a) 2 IN 1 ) 



m 

m 2 + f 2 ’ 


(Bll) 


where T = —1/(2 H p ) + g/c 2 is the anisothermal form 
of Eckart’s coefficient. Since c 2 = ygH, then eliminating 
r between (B5) and (B6) yields 

4> = N 2 yH[a(m, <w) + ub(m, <o)]£, (B12) 


where y is the ratio of specific heats. This replicates 
(23-6) of Gossard and Hooke (1975) and can be reex- 
pressed in terms of the relative pressure perturbation 




y[a(m, a>) + i b(m, w)]^. 


(B13) 


APPENDIX B 

Acoustic-Gravity Wave Polarization Relations 

On splitting atmospheric field variables into the sum 
of a mean and perturbation term (X = X + X'), isolating 
and linearizing perturbation terms, and assuming wave 
solutions of the form 

X' = X exp[t(A:x + Mz — Hr)] (Bl) 

= X exp[i(/:x + mz ~ kit) + z/2H p ], (B2) 


A linearized perturbation expansion of the ideal gas law 
yields 


£ = £ _ £ 

T p p’ 

which, on inserting (B9) and (B13), yields 


(B14) 


L 

T 



(1 _ (y — l )[a{m, a>) + tb{m, <w)]}£ 


(B15) 
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Equations (B3)-(B7) also yield the polarization relation 
-(co 2 - f 2 ) 

= t( w -„■) <" + < B16 > 

For an isothermal ( H = H p ) nonrotating (/ = 0) at- 
mosphere, it can be shown that the dispersion and po- 
larization relations (B8)-(B16) are the same as Eqs. 
(21)-(24) and (26) of Hines (1960). 

On linearizing the tracer continuity equation (lb), as- 
suming a height-varying mean tracer profile n(z), and 
substituting the wave solution (Bl), we obtain 

„ dn _ - 

-turn H w + n V • U = 0. (B17) 

dz 

From (B16), it follows that 

V • U = iku + iMw, (B18) 

1 N 2 1 

~ ——[a{m, (6) + ib(m, <u)] >w, (B19) 

whereupon (B17) yields 



it 

n 


[ 1 1 dn N 2 

Uh + *Tz~ J [a(m ’ a) + Lb{m ’ M)] 


l 


(B20) 


Thus, on comparing (B9), (B13), (B15), and (B19) 
with the simplified parcel-based derivations in appendix 
A (see also sections 2 and 3), we see that the normalized 
terms I a(m, o>)l and I b(m, oj)I must be *^1 for the two 
sets of equations to be equal to a good approximation. 
We note from (BIO) and (Bll) that the frequency de- 
pendence of a(m , co) and b(m , co) is weak, so for sim- 
plicity we consider the limiting maximum values 


tfmaxOO 


c 2 \m 2 + r 


(B21) 


*(«) = 4 


+ r 2 


(B22) 


that arise for ( co/N ) 2 0. 

Since these maximum values depend only on the ver- 
tical wavenumber m, Fig. Bla plots tf max (m) and fc max (m) 
versus A z = 27r/lml for an isothermal atmosphere of T 
= 230 K. Figures Bib and Blc show the corresponding 
maximum difference in magnitude and phase between 
the polarization relations (B9) and (B15) andjthe cor- 
responding parcel approximations p/p = —TIT = (N 2 / 
g)£ given by (18). Other plots (not shown) revealed that 
these differences increase (decrease) with decreasing 
(increasing) T, and that the equality in magnitude of the 
density and temperature departures in Fig. Bib occurs 
only when the atmosphere is isothermal. For vertical 
wavelengths A z ^ 20 km, we see from Fig. B 1 that the 
approximate parcel-based formulas of appendix A, ap- 
proximated by the wave polarization relations (14) and 
(18), produce only small deviations in amplitude and 


phase from the full acoustic-gravity wave polarizations 
relations (B9) and (B15). 


appendix c 


Photochemical Reponse of Ozone Mixing Ratio in 
a Vertically Displaced Air Parcel 

We start with the parcel-based ozone rate equation 
(Lindzen and Goody 1965) 


^ _ 2 j n _ 

I. J 

dt ■ ■ k 0i n 02 

where 


(Cl) 



e(1.3 X 10 28 m~ 3 ) 


300 K 


2.3 


X exp[— (2060 K )/T] = G(T) (C2) 


(Demore et al. 1994). Equation (Cl) differs from (2) 
and (31) only through insertion of an extra factor, e, 
that scales the ozone destruction rate due to recombi- 
nation with monatomic oxygen, k 13 , to a net rate that 
resembles the one obtained when reactions with NO,. 
HO x , and CIO, species are included (Zhu and Holton 
1986). 

We use (Cl) and (C2) to describe the time evolution 
of ozone in a parcel displaced adiabatically by a small 
vertical distance £. Since (A18a) holds for n 0 ,(z 2 )/ 
n o 2 {z 2 ) ~ 1, then 


dq o,fe) 

dt 


2/0,^03(22) 


- -- T* zG(T(z 2 ))qo 3 (z 2 ) exp (~£/H D ), 
n oA<>2) 

(C3) 

where q Q} (z 2 ) = ^ 0 ,(z 2 ) + q' 0i (z 2 ) from (A19). We sim- 
plify (C3) by first noting that the background produc- 
tion/loss rate 


/?o 3 (Z2) = 2/0,70,(22) - 27, 


G(T (z 2 ))q p,(z 2 ) 


o 3 


0. 


(C4) 


We then use a TMS expansion of exp(— £/H D ) and a 
truncated Taylor series expansion 


G(T + T') = G(T ) + T' 



G + G', (C5) 


where T'(z 2 ) = T(z 2 ) — T(z 2 ). Finally, on assuming 
constant q Ql (no diffusive separation), constant J x , and 
retaining only linear terms in the perturbation expansion 
of variables in (C3), we obtain, with the aid of (A 13), 
(C4), and (C5), 
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dq 0 Jz 2 ) 

dt 


-B(z 2 )q' 0 ,(z 2 ) ~ C (z 2 )T'(z 2 ), (C6) 


where 


B(z 2 ) 


4J 0 ^(z 2 )G(T(z 2 )) 

«o 2 fe) 


(C7a) 


Cfe) = 


B(Z2)q(z 2 ) 

2T*(z 2 ) 


T(z 2 ) - 


f d(\nG) \ 


(C7b) 


Equation (C6) then yields the well-known ozone per- 
turbation result of Lindzen and Goody (1965) (see sec- 
tion 3f), with the expressions (C7a) and (C7b) essen- 
tially the same as those derived by Zhu and Holton 
(1986). 1 More complete expressions for R Q3 with the 
NO*, HO*, and CIO* rate terms included explicitly also 
leads to a linearized ozone deviation equation of the 
same form_as (C6), although the formulas for the co- 
efficients B and C are more complicated in this case 
(see, e.g., Stolarski and Douglass 1985). 
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Abstract. Power spectral densities (PSDs) of 
mesoscale fluctuations of temperature and rate of 
change of temperature (heating/cooling rate) due to a 
spectrum of stratospheric gravity waves are derived us- 
ing canonical spectral forms based on observations and 
linear gravity wave theory. The parameterization de- 
veloped here assumes a continuous distribution of hor- 
izontal wave phase speeds, as opposed to a previous 
spectral parameterization in which all waves were as- 
signed stationary ground-based phase speeds. Signif- 
icantly different heating/cooling rate PSDs result in 
each case. The differences are largest at small hori- 
zontal scales, where the continuous phase-speed param- 
eterization yields heating/cooling rate PSDs that are 
several orders of magnitude smaller than in the sta- 
tionary phase-speed parameterization. A simple Monte 
Carlo method is used to synthesize randomly-phased 
temperature perturbation time series within tagged air 
parcels using either spectral parameterization. These 
time series are incorporated into a nonequilibrium, mi- 
crophysical trajectory-box model to assess the micro- 
physical consequences of each parameterization. Col- 
lated results yield a “natural” geophysical scatter of 
instantaneous aerosol volumes within air parcels away 
from equilibrium conditions. The amount of scatter was 
much smaller when the continuous phase-speed param- 
eterization was used. 
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1. Introduction 


The importance of polar stratospheric clouds (PSCs) in ozone depletion has prompted 
considerable research into their composition and formation. Observations show that PSCs 
often form preferentially over mountainous terrain, where they seem to be associated with 
the penetration of mountain lee waves into the stratosphere (e.g. Cariolle et al 1989; 
Deshler et al 1994; Godin et al 1994; Fromm et al. 1997; Carslaw et al. 1998a, b). 
Recent modeling work has shown that mountain wave temperature perturbations affect the 
formation and growth of aerosol particles, as well as subsequent heterogeneous chemical 
reactions which take place on the particles (e.g., Peter et al 1994; Meilinger et al 1995; 
Bormann et al 1997; Tsias et al 1997). This occurs because microphysical processes have 
a strongly nonlinear dependence on both temperature (e.g., Tabazadeh et al 1994; Carslaw 
et al 1994) and the rate of change of temperature (e.g., Meilinger et al 1995; Tsias et 
al 1997). Thus, details in air parcel temperature histories can have an important impact 
on microphysical end products, such as particle size distributions, total aerosol volume and 
particle freezing (Drdla and Turco 1991; Jensen and Toon, 1994; Koop et al 1995, 1997; 
Tsias et al 1997). Furthermore, several important heterogeneous reaction rates have a 
strong temperature dependence (e.g., Bormann et al 1997), which leads to rates of chlorine 
activation in stratospheric mountain waves that are orders of magnitude faster than those 
encountered synoptically (Carslaw et al 1998a). Carslaw et al (1998b) have proposed that 
the combined effect of all such mountain waves in the Arctic could cause an overall decrease 
in ozone levels. The temperature oscillations of a gravity wave have also been shown to be 
important in microphysical models of polar mesospheric clouds (Turco et al 1982; Jensen 
and Thomas 1994) and of high cirrus clouds (Potter and Holton 1995; Jensen et al 1996; 
Lin et al 1998). 

All of these model studies considered the microphysical effects of single waves, whereas 
observations reveal that, on average, broad pseudo-random spectra of temperature fluctua- 
tions exist in the stratosphere (e.g., Nastrom and Gage 1985; Tsuda et al 1991; Bacmeister 
et al 1996; Whiteway et al 1997). Tabazadeh et al (1996) parameterized this by assuming 
Gaussian probability density functions (PDFs) of mesoscale temperature fluctuations due to 
a spectrum of waves, basing their variances on stratospheric data taken from the Microwave 
Temperature Profiler (MTP) on NASA’s ER-2 research aircraft. Murphy and Gary (1995) 
(hereafter MG95) also addressed this issue by deriving mesoscale horizontal wavenumber 
spectra of temperature, heating/ cooling rate (i.e., rate of change of temperature) and PSC 
supersaturation fluctuations, based on the temperature spectra observed from MTP and 
other in-situ aircraft measurements in the lower stratosphere. 

The calculations of MG95 were based on the implicit assumption that mesoscale tem- 
perature fluctuations observed by the aircraft were stationary with respect to the ground, 
implying a spectrum of gravity waves with stationary ground-based horizontal phase speeds 
(such as mountain waves). However, analysis of accumulated Meteorological Measurement 
System (MMS) data from stratospheric ER-2 flights led Bacmeister et al (1996) (hereafter 
B96) to different conclusions. They compared fluctuations in these data with idealized mod- 
els of the two-dimensional (2D) spectrum of gravity- wave-induced horizontal velocities as a 
function of vertical wavenumber m, and intrinsic frequency u>, G UjV {m,uj). This continuous 
2D spectrum implies a broad, continuous distribution of wave phase speeds, in contrast to 
the MG95 model. B96 used linear gravity wave theory to transform G UjV (m,uj) into cor- 
responding spectral predictions for the one- dimensional horizontal wavenumber ( k ) spectra 
of horizontal velocities, F U)V (k), and vertical velocities, F w (k ) (Gardner et al 1993a, b). 
Excellent agreement was found between these theoretical spectra and corresponding spectra 
derived from MMS data, leading them to conclude that mesoscale temperature and velocity 
fluctuations in stratospheric ER-2 data were consistent with a spectrum of gravity waves 
with a broad distribution of horizontal phase speeds. 

In this paper, we go on to develop a spectral parameterization of mesoscale tempera- 
ture and heating/ cooling rate perturbations based on the continuous phase speed model of 
Gardner et al (1993a, b) and B96, and compare it with the stationary phase speed param- 
eterization of MG95. We demonstrate that, in each case, very different temperature and 
heating/cooling rate perturbations arise within air parcels as they are advected through the 
wave fields. Reanalysis of MMS data reveals that PDFs of mesoscale temperature ampli- 
tudes in the stratosphere are approximately Gaussian, with depatures only in the “wings” 
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of the distribution. We associate the Gaussian components of the PDFs with the basic 
background spectrum of wave activity that is always present in the stratosphere, and the 
nonGaussian tails with intermittent large-amplitude wave events. We go on to generate 
synthetic time series of parcel-based temperature perturbations using a simple Monte Carlo 
method, which yields Gaussian temperature PDFs that can be governed by either the MG95 
or B96-based spectral parameterizations. Finally, we use these artificial time series within a 
microphysical box- trajectory model for the stratosphere to assess the microphysical impacts 
of each parameterization. 

2. Spectral Parameterization 

In their analysis of velocity and temperature fluctuation spectra from MMS data from 
the lower stratosphere, B96 used linear polarization and dispersion relations for hydro- 
static, nonrotating gravity waves to transform models for the 2D gravity- wave spectrum, 
G UjV (m,o;), into corresponding predictions for the one dimensional horizontal wavenumber 
spectra of horizontal velocities, F UfV (k), and vertical velocities, F w (k ), that are measured 
by the aircraft (Gardner et al 1993a, b). A further assumption of horizontal isotropy in 
wave propagation directions was made in these theoretical derivations. B96 found excellent 
agreement between the velocity spectra observed by MMS and theoretical predictions for 
F U)V (k) and F w (k) based on a separable analytical hydrostatic gravity- wave spectral model 
of the form 


G u ,v — EqA uv {tti)B uv ( cc >), ( 1 ) 

where A UiV (m ) and B UjV {u) are normalized ID model spectra and Eo is the total wave energy 
per unit mass (Fritts and VanZandt 1993). The shapes of the ID model spectra are discussed 
in depth by Gardner et al (1993a, b), Fritts and VanZandt (1993) and B96. Briefly, the 
vertical wavenumber spectrum A U)V (m ) — ► m ~ 3 at large m and peaks at a characteristic 
wavenumber ra*, which is around 27r (2-3 km) -1 in the lower stratosphere. The normalized 
frequency spectrum takes the form B UiV (u)) oc [1 -f (f / lj) 2 ]cj~ p between the allowed intrinsic 
frequencies \ f \< \ uj \< N (B U)V (lj) = 0 elsewhere), where / is the Coriolis parameter, N is 
the background Brunt-Vaisala frequency, and p = 2. These shapes are based on both theory 
and accumulated data, and although nonseparable models for G UjV (m,u;) have also been 
proposed (e.g., Gardner 1996; Dewan 1997), the separable form (1) is still used commonly 
in the stratosphere (e.g., Gardner et al 1993a, b; Fritts and VanZandt, 1993). The model 
(1) implies a broad continuous distribution of wave phase speeds within the spectrum, as 
opposed to the stationary phase speed model of MG95. 

We now apply the model (1) to predict PSDs of temperature and heating/cooling rate 
perturbations. 

2.1. Temperature PSDs 

B96 also evaluated PSDs of MMS potential temperature data and compared them to the 
predictions of the model (1). Simple polarization relations for hydrostatic gravity waves lead 
to the relation 


G T (m,w)=(^) G@(m,u) = Q(w) (~^-) G UiV (m,u), (2) 

where is the 2D gravity-wave temperature spectrum, Ge(m ,w) is the potential 

temperature spectrum, g is gravitational acceleration, 0 and T are background potential 
and absolute temperatures, respectively, and Q(u>) is the ratio of potential to kinetic energy 
of a given wave of intrinsic frequency to. For hydrostatic waves, Q(u>) = 1 for nonrotating 

(a; 2 — f 2 ) 

flow and Q(u>)= for rotating flow (for the calculations in this paper, we have assumed 

a nonrotating atmosphere). This leads to model horizontal wavenumber PSDs 

F t ( k) = Tq Ft ( k) oc F @ (k) oc (^) F u , v (k), 


( 3 ) 



where T 0 2 is the total temperature variance and F T (k ) is normalized model spectrum. Thus, 
the potential temperature, absolute temperature and horizontal velocity fluctuations should 
all have similarly-shaped horizontal wavenumber spectra. B96 showed that F&(k) and 
F u ,v(k) computed from MMS data had similar shapes, and thus agreed well with the model 
prediction (3) based on (1). 

Of additional relevance to microphysics is the intrinsic frequency PSD of temperatures, 
which governs the time variations of temperature fluctuations within tagged air parcels as 
they are advected by the background flow. This spectrum is difficult to measure directly. 
However, it is easily derived from the model spectra (1) and (2), to yield 

B t {lo) = TqB t (uj) = Q(lo) (J-^j E 0 B UtV (w) = (4) 

where Bt(w) is the normalized temperature spectrum. 

The present analysis can be contrasted to that of MG95, which assumed (implicitly) 
that all the waves in the spectrum were stationary with respect to the ground. MG95 
used a model for the horizontal wavenumber temperature PSD based on data from various 
airborne platforms, which exhibited similar spectral shapes to the model spectrum Ft{F) 
in (3). Under their stationary phase-speed assumption, a Taylor transformation from the 
ground-based to the intrinsic frame 1 yields the intrinsic (parcel-based) frequency PSD 


B$f G9 *(u>,U) 


-,2 F T (u/U) 

0 U 


( 5 ) 


where U — (U 2 -f U 2 ) 1 / 2 is the background wind speed and u = kU. This derivation 
also assumes no directional variability in F T (k) (i.e. azimuthal isotropy). Note the strong 
dependence of this parcel-based spectrum on the background wind speed. 


2.2. Heating/Cooling Rate PSDs 

If the wave spectrum is nondissipating, as is often the case in the lower stratosphere [see, 
e -g*> Figures 7 and 8 of Tsuda et ai (1991)], an air parcel will be displaced adiabatically 
about its equilibrium position by the wave field. Thus, the potential temperature of the 
parcel will be conserved, which, for a single wave, yields the linearized perturbation equation 

D t e f = d t e f + Udj + vdyO' = -Vd,e, (6) 

where 0 and w are wave-induced perturbations of potential temperature and vertical veloc- 
ity, respectively, U and V are the background zonal and meridional winds, respectively, 0 is 
the background potential temperature, subscripted d terms denote local partial derivatives 
(i.e. d t = d/di ), and D t is the intrinsic (parcel-based) time derivative. Horizontal gradients 
in background quantities were neglected in deriving (6). 

Hydrostatic gravity waves have shallow vertical structure, and so the temperature and 
potential temperature perturbations are approximately related, as in (2), according to 

6' V 


Hhe Taylor transformation converts stationary horizontal 
wavelength structures in one frame (in this case, the ground- 
based frame) into temporal variations measured within another 
frame which is being advected through this horizontal variability 
(in this case, the intrinsic frame). Note that this is the inverse 
Taylor transformation to the one used to derive ground-based fre- 
quency spectra due to stratified two-dimensional turbulence mod- 
els, in which the horizontal variability is assumed to be stationary 
(“frozen”) in the intrinsic frame [see, e.g., eq. ( 9 ) of Vincent and 
Eckermann ( 1990 )]. 



(e.g., Eckermann ei al 1998). Thus, rearranging terms in (6) and scaling with (7), we 
obtain the expression 


D t T' = d t T' + Ud x T' + Vd y T' = -w'd z Q (| ) = - (^) , (8) 

for the Lagrangian rate of change of temperature perturbation (i.e., the heating/cooling rate 
perturbation within an air parcel advected by the mean wind). Thus, for linear adiabatic 
gravity waves, the heating/cooling rate perturbations are proportional to the vertical velocity 
perturbations, and so the heating/cooling rate PSD, Fj) i 'r{k) ) follows from simple scaling 
of the PSD of vertical velocity according to (8): 

F DtT (k) = (^) F w (k) = FE$(k). (9) 

The horizontal wavenumber PSD of heating/ cooling rate fluctuations under the continuous 
phase speed model, (9), may be contrasted with the one arising from the MG95 stationary 
phase speed model, which is given by 

F™$ 95 (k,U) = U 2 k 2 F T (k). (10) 

This relation follows via a Taylor transformation from the ground-based frame (where 3* 
terms vanish since waves are taken to be stationary) to the intrinsic frame to derive the 
parcel-based time derivative D t (Murphy and Gary 1995). Note that, unlike (9), the MG95 
heating/ cooling rate PSD (10) is a strong function of the background wind speed U. 

Intrinsic frequency PSDs of heating/cooling rate follow from the corresponding tempera- 
ture PSDs. For the B96-based model, 

Bd^t( u ) — u) 2 B^ 96 (lo) } ( 11 ) 

and so is approximately flat given the approximately u>~ 2 shape of B^ 96 (lo) implied by (4). 
The corresponding PSD from the MG95 model is 

B%% 95 ( u,U) = u 2 b¥' G95 (u,U) = j^F t (u/U), (12) 

where to = kU. 


3. Comparisons with Data 

To compare the various parameterizations with data, we use vertical velocity and temper- 
ature fluctuations derived from MMS data during 3 aircraft campaigns employing NASA’s 
ER-2 aircraft: (1) the second Arctic Airborne Stratospheric Expedition (AASEII); (2) the 
Stratospheric Photochemistry Aerosol and Dynamics Experiment (SPADE); and (3) the 
Airborne Southern Hemisphere Ozone Experiment /Measurements for Assessing the Effects 
of Stratospheric Aircraft (ASHOE/MAESA) campaign. Locations and dates of these cam- 
paigns, as well as a brief description of the MMS data set, can be found in B96 and refer- 
ences therein. The PSDs used in this study are calculated from discreet wavelet transforms 
(DWTs) in the manner outlined in B96. 

3.1. Heating/Cooling Rates PSDs 

It was noted in section 2 that the two parameterizations use a similar observationally- 
constrained model for the horizontal wavenumber spectrum of temperatures, Er(fc), but 
derive different heating/cooling rate PSDs from it. These differences are ilustrated in Fig- 
ure 1. The solid curve with triangular symbols in Figure la shows Fp 9 j,(k) derived from 
MMS vertical velocity spectra using (9). It represents a log-average over all flights in the 
AASE II, SPADE, and ASHOE/MAESA campaigns. The remaining curves in Figure la 
show F^ t ^ 95 (k,U) for three U values of 1 m s -1 , 20 m s _1 and 50 m s -1 , as obtained by 



scaling a corresponding log-averaged PSD of MMS temperature (shown in Figure lb) using 

The shapes of the two heating/cooling rate PSDs are clearly very different. The MG95 
PSDs are much flatter at small scales. At scales smaller than 6.4 km, F$1}*\k y U) decays 
weakly with a spectral index between 0 and —1, while Fp^^k) decays more rapidly with a 
spectral index close to -3. The magnitude of F^ i ^{k ) U) also depends sensitively on the 
background flow speed U , since it determines how rapidly the parcels are advected through, 
and perturbed by, the stationary horizontal wave structures. In contrast, the heating/cooling 
rate PSD determined from the B96-based parameterization (9) has no explicit dependence 
on background flow conditions. Observationally, B96 found a weak dependence of the MMS 
vertical velocity PSD on the magnitude of the background horizontal wind, and other studies 
have also reported some sensitivities of stratospheric velocity spectra to the background 
wind profile (e.g., Eckermann 1995). The sensitivity of the MMS vertical velocity PSDs 
to background wind speed appeared to vary with horizontal scale: at k>m* the vertical 
velocity PSDs appeared to be less sensitive to changes in the background wind speed than 
the PSDs at k<m * (Bacmeister et al 1996). So, while some evidence for a dependence on 
wind speed exists, existing data indicate that it is significantly weaker than the quadratic 
dependence of F^ 95 (k,U) on U in (10) and Figure la. 

Table 1 gives spectrally-integrated root-mean-square (rms) heating/cooling rates for the 
B96-based heating/ cooling rate spectrum within several logarithmically-spaced wavenumber 
ranges, as determined by numerically integrating F%$(k) in Figure la. Integrated heat- 
ing/cooling rates based on (9) peak at scales between ~3-25 km. This is in contrast to the 
MG95 spectral model, which yields integrated cooling rates that increase monotonically as 
horizontal scales decrease (see Table 1 of MG95). 

3.2. Temperature PDFs 

Another approach to parameterizing multi-wave temperature perturbations within micro- 
physical models is to employ Gaussian PDFs (Tabazadeh et al 1996). The DWT analysis 
of MMS data by B96 preserves information about the location and amplitude of fluctuations 
at each of the characteristic scales (or “octaves”) resolved by the DWT, which allows us to 
study temperature fluctuation statistics and PDF’s as a function of horizontal scale. 

Figure 2 shows PDFs of DWT temperature amplitudes at different horizontal scales. 
Least-squares fits to a Gaussian distribution over the central 68% of the data are also 
shown. Within ±la of the peak, the PDFs closely approximate a Gaussian distribution, as 
assumed by Tabazadeh et al (1996). The rms amplitudes of the fluctuations in Figure 2 
are comparable to those obtained by MG95 at similar horizontal scales. 

However, there is evidence of non-Gaussian high-amplitude tails in all of the PDFs in 
Figure 2. In other words, high-amplitude temperature fluctuations are somewhat more 
probable than would be inferred from the Gaussian distributions determined from fits to the 
central 68% of the data. These departures from Gaussianity become more pronounced as the 
horizontal scale of the fluctuations decreases in Figure 2. Chi-squared tests of the observed 
PDFs versus the Gaussian fits indicate that the non-Gaussian behavior in the wings of the 
distributions is highly significant statistically. Similar behavior is observed in the PDFs of 
horizontal velocity fluctuations measured by MMS in the stratosphere (Bacmeister et al 
1997a). The origin of the non-Gaussian tails in the temperature fluctuations is not entirely 
clear. However, intermittent large-amplitude wave events are likely to be at least partially 
responsible (Bacmeister et al 1990; Alexander and Pfister 1995). 

4. Parameterizations for Parcel-based Calculations 

The spectral parameterizations (3), (9), (10) and those of MG95 are presented as hori- 
zontal wavenumber spectra, which govern the horizontal “line” distribution of a collection 
of parcels at a given instant in time. However, many microphysical models are Lagrangian, 
concentrating instead on time variations of meteorological and microphysical parameters 
within tagged air parcels as they are advected by the background flow. 

For such models, spectral parameterizations of wave-induced temperature fluctuations 
following an air parcel are needed. These are given by the intrinsic frequency spectra (4) 



and (5) for the B96-based and MG95 parameterizations, respectively. Again, the two pa- 
rameterization give very different results. For temperatures, quasi-invariant spectra ~oj ~ 2 
arise in the B96-based model (4), whereas the MG95-based spectrum (5) follows the shape 
of the horizontal wavenumber spectrum according to Fr(u;/7/), and its intensity scales with 
background wind speed U . Analogous differences exist between the heating/cooling rate 
spectra (11) and (12). 

To assess the effects of each parameterization on Lagrangian microphysical models, it is 
first necessary to generate parcel-based temperature perturbation time series governed by 
these spectral models. A straightforward approach for the B96-based parameterization is to 
generate a random time series of the form 

^£96 00 — V^To y ' sin(u;^ -f <^ 2 -), (13) 

2 — 1 

where y/2 transforms from rms to peak amplitudes, Au; is a (constant) frequency step, 
uji is the ith. (evenly-spaced) frequency and <j)i is a random phase for this harmonic. The 
frequency step Au; must be small enough to resolve shape transitions in By 96 (u>) adequately, 

and n and Au; are chosen so as to integrate B^ 96 (ui) over its full range between / and N 
(i.e., uq = / + Au;/2, w n = N — Acj/ 2). While this is similar to the Monte Carlo method of 
Eckermann (1990), many other randomly-phased Fourier and/or wavelet synthesis methods 
could also be used (e.g., Owens 1978; Benzi et al 1993; Alexander 1996; Elliott et al 
1997). Note that (13) yields Gaussian temperature PDFs, and so is also consistent with the 
PDF-based modeling strategy of Tabazadeh et al (1996). 

The main remaining task is to choose a representative total rms temperature amplitude To 
for insertion into (13). Here, we attempt a data-based estimate of it from our temperature 
PSD in Figure lb, which was computed over a limited wavenumber range &i<fc<& 2 , where 
ki = 27r(102.4 km) -1 and k^ = 27t( 0.8 km) -1 . Thus, integrating this spectrum gives us some 
fraction, 6, of Tq . We estimate e by numerically integrating Gt(^,u;) in (2) over all m and 
uj and comparing it with a second integral constrained within the range k\<k<k2. Figure 3 
shows a sample contour plot of Gr(ra,u;) for m*= 27r(3 km) -1 . The shaded region shows 
the integration range on restriction to the horizontal wavenumber range Using 

a typical stratospheric value of ra*= 27r(3 km) -1 (Tsuda et al 1991), we obtain e = 0.21 
using representative midlatitude / and N values. Thus, according to the continuous phase 
speed model (2), a little over one fifth of the total variance T 2 is contained in the horizontal 
wavenumber PSDs shown in Figure lb and B96. 

On numerically integrating the ID temperature PSD in Figure lb, we obtain a variance 
of 0.16 K 2 , and thus a total variance T 2 = 0.16/e «0.8 K 2 . This is a small value compared 
to other measurements and analyses (e.g., Tsuda et al 1991; Allen and Vincent 1995), but 
it should kept in mind that a significant fraction of the total ER-2 flight miles during these 
experiments occurred over open ocean, where wave variances are smaller than those over land 
(e.g., Nastrom et al 1987). Furthermore, the mean spectra in Figure lb and B96 are log- 
averages of wavelet spectra taken during 1024 s flight intervals. Log-averages deemphasize 
periods of high activity, and thus the variance of the log-average spectrum is less than the 
mean variance of the data. Nevertheless, our estimate of Tq« 1 K 2 can be considered broadly 
representative of “ambient” stratospheric conditions, or a “lower bound” estimate. In any 
given situation, Tq may be considerably larger, and may depend on environmental factors 
such as underlying cloudiness and terrain, background winds (at all levels), altitude, and so 
on. A sample parcel-based temperature fluctuation time series using T 2 = 0.76 K 2 in (13) is 
shown in Figure 4a. 

Temperature perturbation time series for the MG95 parameterization can also be gen- 
erated by replacing B^ 96 (oji) in (13) with the corresponding normalized MG95 spectrum, 
W -1 JPr(w</W), from (5). Figure 4b shows the resulting time series, using the same n, To and 
(j)i values as in Figure 4a, U— 50 m s -1 , and where the shape and wavenumber range of the 
temperature PSD in Figure lb was used to specify F T {ui/U). Note the smaller time range 
and the “rougher” appearance of this time series, due to higher frequency content and a 
different spectral shape. This time series also highlights a serious shortcoming of the MG95 
parameterization: that the gravity- wave frequency limits | / |<| to |< N are not imposed 



in this parameterization, and thus unphysical high-frequency temperature oscillations are 
simulated within air parcels when U is large. 

An additional general caveat must be made about the use of (13) to generate parcel- 
based temperature fluctuations. For multi-parcel trajectory models, we would expect some 
degree of spatial correlation between the temperature perturbations within closely adjacent 
parcels. A straightforward application of (13) within such models will not provide these 
correlations. However, (13) can be generalized to provide spatio-temporal parcel-based 
temperature fluctuation fields T (x, y, z, t) using the complete four-dimensional model power 
spectrum ^^(Ar, /, m, a;), assuming isotropy or some other simple directional dependence in 
wave propagation directions (e.g., Gardner et al 1993a 1993b). 

5. Microphysical Effects of the Parameterizations 

Figures 4a-b show that the B96-based and MG95 parameterizations lead to noticeable 
differences in simulated mesoscale temperature variability along an air parcel trajectory. 
Meilinger et al. (1995) and Tsias et al. (1997) have shown that rapid temperature fluctua- 
tions can cause the composition of liquid stratospheric aerosols to depart from equilibrium. 

The extent of non-equilibrium depends upon the cooling and heating rates of air parcels 
(Tsias et al 1997), with more rapid cooling and heating leading to greater departures from 
equilibrium. Under conditions where the stratospheric aerosol is composed principally of 
H2SO4 and H2O (typically at temperatures T > 195-200 K), the effect of temperature fluc- 
tuations can largely be ignored since equilibrium between gas phase H2O and the droplets 
is rapidly achieved. The greatest departures from equilibrium occur under conditions where 
HNO3 is strongly partitioned from the gas phase into the liquid aerosol, which is typicallv 
at T < 195-200 K (Carslaw et al 1997). 

Here, we investigate how each parameterization affects the evolution of ternary H2O/HNO3/H2SO4 
aerosol droplets by inserting the time series in Figure 4a-b into the nonequilibrium micro- 
physical box-trajectory model of Meilinger et al (1995). We also investigate the effect 
of these temperature fluctuations on the growth and evaporation of nitric acid trihydrate 
(NAT) particles. Our aim is to examine how the different natural background temperature 
fluctuations might affect the instantaneous volumes of these typical stratospheric particles. 

In the experiments using this model, parcels were advected along the 0 = 454 K isen- 
trope using an initial water vapor mixing ratio of 5 ppmv and an initial gas-phase HNO3 
mixing ratio of 10 ppbv. The perturbation time series in Figure 4a-b were superimposed 
upon the mean temperature, and four separate experiments were conducted for each time 
series using T= 188, 190, 191.5 and 193 K (with corresponding mean pressures p=47, 48.7, 

50, and 51.4 hP, respectively 2 ). Total aerosol volume and temperature within the parcel 
were calculated and recorded at each time step. Results for all four experiments using the 
B96-based and MG95-based temperature perturbation time series are collated and plotted 
in Figures 4c and 4d, respectively. The thin black line on both panels represents the thermo- 
dynamical equilibrium curve of Carslaw et al (1994). Also indicated on each figure is the 
approximate amount of scatter in observed aerosol volumes (in terms of the equivalent scat- 
ter in temperature) measured by Dye et al (1992) on 24 January 1989. These observations 
were used by Carslaw et al (1994), Tabazadeh et al (1994) and Drdla et al (1994) for 
comparison with volumes calculated using a liquid aerosol model assuming thermodynamic 
equilibrium. 

Figures 4c and 4d show that the temperature fluctuations in both the B96-based and 
MG95 parameterization lead to significant scatter in aerosol volume around the thermody- 
namic equilibrium value. This is because the time constant for diffusive growth and shrinking 
of the droplets by HNO3 mass transfer is rather large compared to the rate of change of air 
parcel temperature that induces the size changes [for a discussion of these time constants, 


2 Note that Carslaw et al. (1994) used a constant pressure 
of 55 hPa in calculating aerosol volumes as a function of tem- 
perature, implying a change in potential temperature between 
the different temperatures. Here, the calculations refer to aerosol 
volumes on a constant potential temperature level. The difference 
in volume is small. 



see Meilinger et al (1995)]. In particular, it is the large droplets whose growth is diffusively 
hindered, while, as pointed out by Meilinger et al (1995) and Tsias et al (1997), the small 
droplets adjust their mass rapidly to changes in temperature. However, these small droplets 
do not contribute significantly to the total aerosol mass. For both the B96-based and MG95 
gravity- wave parameterizations, the departures in aerosol volume from equilibrium are even 
greater than those produced by rapid synoptic-scale cooling (Carslaw et al 1997). 

The results using the MG95 parameterization (Figure 4d) show a great deal more scatter 
about the equilibrium line in the 189 to 193 K range than the B96-based results in Figure 4c. 
Since there is more rapid heating and cooling of parcels in the MG95 parameterization (see 
Figure 4a-b), large aerosol droplets can be rapidly heated to high temperatures before they 
can adjust by evaporating mass, leading to larger departures from the equilibrium curve. 
Aerosol volumes a factor of 5-10 larger than equilibrium predictions occur in Figure 4d. 
Nonequilibrium behavior is also evident in the B96-based experiments (Figure 4c) but is 
much less pronounced, with aerosol volumes generally less than a factor of two greater than 
equilibrium predictions. The strong nonequilibrium behavior in the MG95 parameterization 
also produces a higher total mean aerosol volume than in the B96-based parameterization. 

Clearly the amount of scatter in the B96 case is much closer to that observed by Dye et 
al (1992) than is the case for the MG95 trajectory. This may reflect the fact that the MG95 
trajectory yields unrealistically rapid cooling and heating rates, and therefore unrealistically 
broad scatter of calculated aerosol volumes. The calculated scatter in liquid aerosol volumes 
in the B96 case indicates the “natural” range of scatter along an individual trajectory arising 
from a background of mesoscale temperature fluctuations. The scatter observed by Dye et 
al (1992) most likely includes this effect, but also contains scatter arising from air-mass 
variations, as well as from instrument noise, particularly at low volumes. Another source of 
scatter in the measurements is that due to different particle types (most likely in different air 
masses). The Dye et al observations from 24 January 1989 appear to be consistent with pure 
liquid ternary aerosols at all temperatures, which is not the case for observations on other 
days. Note that the calculated scatter shown here is apparent only where aerosol volumes 
increase sharply with decreasing temperature, while scatter is apparent at all temperatures 
in the Dye et al data. From their measurements alone, it is not possible to say how much of 
the observed scatter in aerosol volumes was produced by mesoscale temperature fluctuations. 
However, the observed scatter in the data should represent an upper limit on estimates of 
the true value of scatter due to temperature fluctuations. 

We now consider the effect of the B96-based temperature fluctuations on NAT growth 
(Fig. 5) using an extended version of the Meilinger et al (1995) model (Tsias et al 1997). 
We show a case where 1% of the background liquid aerosols are assumed to be nucleated as 
NAT (corresponding to about 0.1 cm”*). NAT particles were assumed to form and evaporate 
at the NAT equilibrium temperature. It is likely that NAT formation may require a large 
degree of supercooling, or may even occur heterogeneously on ice at temperatures several 
Kelvins lower. The calculated scatter shown here is therefore a lower limit to that expected 
in reality. Figure 5 shows the scatter in total aerosol volume, and is therefore comparable 
with aerosol volume measurements. The effect of liquid aerosols is relatively minor, so the 
volumes shown are also rather similar to those of NAT alone. The amount of scatter is now 
considerably greater than was the case for liquid aerosol growth. The reason is that, with 
only 0.1 NAT particles cm -3 , the individual particles grow rather large and are therefore 
more strongly diffusively-limited than the many smaller liquid aerosols. The amount of 
scatter is the same as the amplitude of the temperature fluctuations themselves. With 100% 
NAT activation, the scatter again resembles that for liquid aerosols, while for even fewer 
NAT particles it is somewhat greater. For the MG95 temperature fluctuations and 1% NAT, 
the amount of scatter is somewhat larger than that obtained with the B96-based time series, 
although the difference is not as pronounced as it was for the liquid aerosols. As the number 
density of NAT particles increases, the difference in scatter between B96 and MG95 NAT 
aerosol volumes also increases. 

These calculations indicate that, even under “passive” background synoptic conditions, 
a large amount of scatter in aerosol volume can occur due to mesoscale temperature per- 
turbations. Although we have compared the calculated scatter with a single set of in situ 
particle volume observations, we believe such a “natural” scatter will be present in all such 
particle ensembles. This sets a natural limit on the “quality” of particle volume and com- 
position measurements. Under some conditions this may make identification of NAT and 



liquid aerosols in the stratosphere from volume observations alone very difficult unless non- 
equilibrium growth is taken into account. 

6. Summary and Discussion 

The work of MG95 indicated that the temperature perturbations produced by a spectrum 
of stratospheric gravity waves had important effects on stratospheric microphysics. Here, 
we presented revised models of the PSDs of mesoscale temperature and heating/cooling 
rate perturbations due to such waves, based on the spectral data of B96. The shapes 
of these PSDs were significantly different from those obtained by MG95. We derived a 
heating/ cooling rate PSD that is more steeply decreasing at the smallest horizontal scales, 
and is far less sensitive to variations in background wind speed (Figure la). Overall, RMS 
cooling rates using the B96-based model were significantly smaller than those derived by 
MG95. 

We developed a simple method for generating random time series of temperature pertur- 
bations within air parcels, as governed by either the MG95 or B96-based spectral param- 
eterization. These time series were imported into a microphysical box model to assess the 
microphysical consequences of each parameterization. Both parameterizations generated to- 
tal aerosol volumes as a function of air-parcel temperature that departed significantly from 
equilibrium (Figure 4c-d). The MG95 parameterization generated aerosol volumes that ex- 
ceeded equilibrium values by up to a factor of 10. The B96-based approach yielded smaller 
departures, with aerosol volumes generally within a factor of 2 of equilibrium predictions. 

The ER-2 aerosol volume versus temperature data of Dye et al (1992) show significant 
scatter about the equilibrium curve. In flights for which air mass variations were estimated 
to be small, aerosol volume as a function of temperature does not appear to show the large 
departures from equilibrium predicted by the MG95 approach (Figure 4 d). As dicusssed 
in section 4, the total temperature variance of T 0 2 = 0.76 K 2 used in Figure 4 is probably 
a lower bound, and use of a larger T 2 in Figure 4 would yield even greater departures 
from equilibrium using the MG95 parameterization. Thus, the degree of scatter in the 
aerosol volume data of Dye et al. (1992) is more closely approximated by a microphysical 
model using a B96-based parameterization of typical stratospheric mesoscale temperature 
variability (Figure 4c). 

The differences between the MG95 and B96-based approaches to parameterizing adiabatic 
gravity-wave heating and cooling perturbations arise from the frequency spectrum assumed 
for the gravity waves in each approach. MG95 assumed (implicitly) a wave field composed 
entirely of stationary waves, which yields a single intrinsic phase speed that is of similar 
magnitude to the background stratospheric wind. The B96-based approach assumes that a 
broad distribution of intrinsic wave phase speeds is present, with a peak centered ~ 10 ms " 1 
in the lower stratosphere. Thus, waves in B96 are assumed to ’’run with” the background 
flow to a large degree, reducing the relative speed of air parcels with respect to the waves 
in regions where background winds are large (e.g., near the edge of the polar vortex). 

The full implications of these differences for stratospheric microphysics are not yet clear. 
Generally speaking, however, both instantaneous and time-averaged aerosol volumes and 
surface areas will be larger using the MG95 parameterization. However, the duration of 
cold temperature events due to gravity waves are longer using the B96-based parameteriza- 
tion. This extended residence time in cold conditions could, for example, lead to increased 
settling of particles. The time evolution of heterogeneous reactions with strong tempera- 
ture dependences could also depend on the shape of the frequency spectrum assumed for 
wave-induced temperature fluctuations. More work is needed to assess these possibilities 
fuHy. 

We consider it unlikely that the background temperature fluctuations described here could 
induce nucleation of recognized nitric acid hydrates. The work of Tsias et al. (1997) shows 
that large amplitude (>6 K peak to peak) rapid cooling and warming events, such as found 
in localised mountain-induced gravity waves, would be necessary. In such events, the compo- 
sition of HNO 3 /H 2 SO 4 /H 2 O droplets can approach almost binary HNO 3 /H 2 O with HNO 3 
concentrations in the liquid as high as 59 wt%, leading to the possible nucleation of nitric 
acid dihydrate from the droplets. Temperature fluctuations as large and rapid as these 
are not found in the gravity wave temperature time series considered here (see Figure 4a). 



Therefore, hydrate nucleation in temperature fluctuations may well be restricted to regions 
of the stratosphere directly affected by large- amplitude mountain waves [see, e.g., Fig. 2 of 
Carslaw et al (1998b)]. Note, however, that these arguments refer only to homogeneous 
freezing of nitric acid hydrates from the liquid. As shown by Tabazadeh et al (1996), 
background gravity wave temperature fluctuations may sometimes be sufficient to push air 
parcel temperatures sufficiently below the ice frost point to cause ice nucleation. However, 
the magnitude of the associated heating/cooling rate fluctuations ( D t T f ) probably have little 
influence on this mechanism. 

From a purely physical viewpoint, we believe that the B96-based parameterization will be 
more realistic than the MG95 approach in most circumstances. In “passive” environments 
away from strong wave sources, gravity wave spectra have characteristic shapes. While the 
processes responsible for generating these spectral shapes are still being debated (e.g., Hines 
1991; Gardner 1996; Dewan 1997), canonical spectral models that are constrained by the 
gravity- wave dispersion relation, like (1), do a good job in interrelating a wide range of 
observed gravity wave spectra. In particular, B96 showed that the shapes and intensities 
of stratospheric horizontal wavenumber spectra measured by aircraft were well predicted 
by these models. The B96-based parameterization presented here is an extension of these 
same models to temperature and heating/cooling rate spectra. However, if spectral models 
other than (1) are preferred, corresponding temperature and heating/cooling rate PSDs can 
be easily derived from them using similar derivations to those in section 2. The MG95 
parameterization, on the other hand, is not constrained by the gravity wave dispersion 
relation, and can generate unphysical high frequency (u> > N ) temperature fluctuations 
within parcels advected by strong winds. 

In regions close to strong wave sources, it is likely that neither spectrally-based approach 
will be the best way to estimate wave temperature perturbations. Alternative models of 
direct wave propagation from orographic or convective source regions show more promise in 
this regard (e.g. Bacmeister et al 1994, 1997b; Alexander 1996; Leutbecher and Volkert 
1996), and have started to be incorporated into regional microphysical models (Carslaw et al 
1998a, 1998b, 1998c). A complete parameterization of gravity wave effects for chemical and 
microphysical models may require both a spectral approach, to capture the small- amplitude 
relatively-ubiquitous background wave field, and explicit event-by-event modeling to capture 
the rarer large- amplitude wave events generated by strong underlying sources. 
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RMS heating/ cooling rates within logarithmically-spaced horizon- 
tal scale ranges, using the B96-based model 
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Figure 1 . (a) Power spectral density (PSD) of heat- 
ing/ cooling rate as a function of horizontal wavenum- 
ber k. (b) PSD of temperature as a function of k . 
PSD’s in (a) and (b) are obtained from log-averages over 
all flights in AASE II, SPADE, and ASHOE/MAESA. 
Filled triangles in (a) show cooling rates obtained by 
scaling the mean PSD of vertical velocity as shown in 
(3). Solid, dotted and dashed lines show cooling rates 
obtained from the horizontal derivative of the mean 
temperature in (b), following MG95. Results shown 
are for 3 different values of mean wind U\ 20, 50 and 
1 m s -1 . 



Figure 2. Logio of normalized PDF’s of tempera- 
ture fluctuation amplitude determined using discreet 
wavelet transforms of MMS temperatures measured 
during ASHOE/MAESA. Thick dashed curves show 
Gaussian fits (which appear parabolic here) to the cen- 
tral 68% of data. Results are shown for fluctuations 
with characteristic horizontal scales of (a) 12.8 km, (b) 
6.4 km, (c) 3.2 km, and (d) 1.6 km. Fluctuation am- 
plitudes appear to be well described by a normal dis- 
tribution in the central or small amplitude portion of 
the distribution. However, PDF’s show non-Gaussian, 
high-amplitude tails especially at smaller horizontal 
scales. These may be associated with infrequent high- 
amplitude wave events, such as waves generated by in- 
tense convection or flow over mountains. 


Figure 3. Model power spectral density (PSD) of tem- 
perature as a function of vertical wavenumber m and 
intrinsic frequency a;, GT(m,u>), from (2) based on the 
separable model (1) of Gardner et al. (1993a). Contours 
of this two-dimensional PSD for m*= 27r(3.0 km) -1 
are shown with solid curves. Cusps in the contours 
at m—m * mark the shape change in A uv {m) in (1). 
Gt(?^,u;) is identically zero for Cc?> | A/" |=0.02 rad s -1 
and u;<|/|=10“ 4 rad s -1 . Shaded region shows por- 
tions of the (m,cj)-plane enclosed by curves of constant 
horizontal wavenumber corresponding to the range of 
horizontal scales between 102.4 km and 0.8 km. For 
ra*=27r(3.0 km) -1 approximately 21% of total variance 
lies within this shaded region. 


Figure 4. (a) Temperature as a function of time for an 
air parcel advected through a random field of gravity 
waves with an intrinsic frequency spectrum B{uj)^uj~ 2 
between / and N (B96-based approach), (b) Temper- 
ature history for an air parcel traversing a stationary 
wave field with a temperature power spectrum like that 
in Figure lb (MG95 approach). A relative speed of 50 m 
s -1 is assumed for the air parcel with respect to wave 
field. A total temperature variance T 2 «0.76 K 2 was 
used in both cases. Total aerosol volume as a function of 
temperature for the two different parameterizations of 
gravity wave temperature fluctuations, (a) and (b), are 
shown in panels (c) and (d), respectively. Aerosol vol- 
umes were calculated for ternary HNO3/H2SO4/H2O 
droplets using the box model of Meilinger et al. (1995). 
Each of the aerosol volume plots in (c) and (d) combines 
the results of 4 calculations using the temperature per- 
turbation histories in (a) and (b), with different mean 
temperatures of 188, 190, 191.5 und 193 K. The thin 
black line in each panel is the thermodynamic equilib- 
rium result for ternary droplets (Carslaw et al. 1994). 
Arrows depict the approximate range of scatter in the 
data of Dye et al. (1992). 



Figure 5. Total aerosol volume as a function of tem- 
perature for NAT particles, where 1% of the background 
liquid aerosols is assumed to be nucleated as NAT. Re- 
sults were derived using the parcel-based temperature 
time series in Figure 4a (B96), using four different mean 
temperatures of 188, 190, 191.5 and 193 K. Solid line 
is NAT volume at equilibrium, dashed line is liquid 
HNO3/H2SO4/H2O droplet volume at equilibrium. 
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Abstract. A two-dimensional model was used to advect 
air parcels under the influence of a linear superposition of 
hydrostatic gravity waves. As wave amplitudes increased, 
vertically profiled wave perturbations became nonsinusoidal, 
producing increased spectral power at large wavenumbers 
M. At saturation amplitudes just before onset of convective 
instabilities, these “tail spectra” assumed universal M~ 3 
shapes, similar to observed spectra in the atmosphere and 
oceans. 

Introduction 

Autospectra of gravity-wave-induced horizontal velocity 
and relative temperature and density fluctuations have char- 
acteristic shapes in the atmosphere and oceans. In par- 
ticular, vertical wavenumber spectra F(M) have a quasi- 
universal M“ 3 form at large M [Smith et a/., 1987]. Many 
theories have been proposed to explain these observations. 
Here, we investigate a specific question of general relevance 
to all such theories: how do the vertically-profiled waveforms 
measured by these sounders relate to the wave-induced ad- 
vection patterns that produced them? 

For a single hydrostatic gravity wave, the answer is 
straightforward: the vertical profile has a sinusoidal form 
consistent with the advection field of the wave. This holds 
even for large-amplitude waves because wave displacements 
remain approximately transverse to the wavenumber vector 
K. However, when a spectrum of J waves superpose, dis- 
placements can no longer be transverse to all wavenumber 
vectors Kj ( j = 1-J) unless they are all parallel. For non- 
infinitessimal amplitudes, it is unclear how this will affect 
measured vertical perturbation profiles. We study this here 
using a simple kinematic model of wave advection. 

Two-Dimensional Advection Model 

A regular spatial grid (x ix1 z iz ) defines the unperturbed 
positions (xf x , z?) of an array of air parcels Pi x ,i z , where 
i x and i z are horizontal and vertical grid indices, such that 

Xi x XQ -f- Zx Ar, i x — 1, ..., Ixf (l) 
Z f z = Zi z = z 0 + i z Az , i* = 1,..., /*, (2) 

and (Ax, A z) is the spatial resolution. We assume a stably- 
stratified inviscid nonrotating background atmosphere of po- 
tential temperature O (z) and no mean flow. 
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We perturb these parcels using a linear superposition of J 
nondissipating hydrostatic gravity waves, which yield parcel 
displacements of the form 


C (Xi x jZi ZJ t') ^ ^ Cj ( Xj x , Zj z , t), (3) 

3 = 1 
J 

- X] & si n (kjX ix + mjZi z - ujt + fa), 


i = i 

j 


X(x ix ,z iz ,t) = ^2 


-rrij 


Cj {.Xi x ? Z iz i £)• 


( 4 ) 


3 = 1 


The subscript j denotes thejth. wave, (j is it’s peak vertical 
displacement amplitude, Kj = (kj , 0, rrij ) is it’s wavenumber 
vector, u)j is it’s oscillation frequency, and is phase. Hor- 
izontal displacements (4) follow from (3) using polarization 
relations [Gossard and Hooke , 1975]. At time t the array of 
parcels Pi x ,i z are advected to perturbed positions 


Xi x + X (Xi x , Zi z , t), (5) 

= Z iz “t" C i.Xi x 5 Zi z , t) , (6) 

yielding an irregularly gridded wave displacement field 

x’ P {xL, iz (t),zf Xyiz (t),t) = X (x ix ,z iz ,t), (7) 

£ ('Ux.tz (i)> Zix,iz (0> 4 = C { X ix i Z iz j t) • (8) 


We resample these displacement fields at the regular grid 
positions (xi x , z% z ) using standard algorithms (Delaunay 
triangulation and interpolation). This yields the array of 
parcels Qi x ,i z (t) that have been advected to these grid po- 
sitions at time £, i.e., 


— Xix 
(*) = zi * 


= + x , (9) 

= i w.(*) + C {*? x ,i z (t),z? x , iz (t),t) . (10) 


The wave displacement fields at these positions are 


X Q (x ix ,z iz ,t)=x ( X ? x , iz (t),z? x>iz (t),t) , (11) 

C Q (x ix ,z iz ,t) = (' (x? x iz (t),z? x iz (t),t) . (12) 

Thus, on the regular grid (x ix ,z iz ), C ( x ix ,z iz ,t ) in (3) 
defines the Lagrangian vertical displacements of the parcels 
Pix,i z from their equilibrium positions (x ix ,z iz ), whereas 
C { X i x i Zi z , t) in (12) defines the Eulerian vertical displace- 
ment field due to parcels Qi x ,i z (t ) that have been instanta- 
neously advected to positions (x ix ,z iz ) at time t. It is these 
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(b) Parcel Displacements at t = 0 



Figure 1 . (a) Wave-perturbed isentropes. Shaded region 
at z — 13.5—14.5 km is expanded in (b), showing advected 
parcels Pi x ,i z (diamonds) and lines connecting them to their 
undisturbed positions ( Xi x , Zi z ). Vertical arrow depicts the 
“beam” of a ground-based sounder at x = 0. 


latter parcels that are sampled by a ground-based profiler: 
e.g., a vertical profiler at x = x ix will measure the vertical 
displacement profile 


Z ( z >t) = C Q (Xix,Z,t). (13) 

The waves also advect potential temperature surfaces. 
This yields a perturbation potential temperature field 


e' Q _(z,t) 

e(z) 


—N 2 (z) 

9 


Z Q (z,t), 


(14) 


where N 2 (z) is the squared background Brunt-Vaisala fre- 
quency [Eckermann et a/., 1998]. Similar fractional relations 
govern the perturbed temperature and density fields. 


Model Results 

Using this model we studied vertical displacement pro- 
files Z (z,to) = £ (xi x ,z,to), as would be sampled at 
some time to by a ground-based vertical profiler at x = x ix . 
We compared Z (z,to) with the corresponding Lagrangian 
perturbation profile Z\z,t) = C* (x ix ,z,t 0 ). Many different 
model simulations were conducted using a variety of choices 
for J, rrij, kj , and <f>j in (3) and (4), but always confin- 
ing \rrij | within some range rrif-mc . The following general 
characteristics were noted. 

The power spectrum of Z* (z,t 0 ), F z (M,t 0 ), was con- 
fined within the wavenumber range M — rrif-mc , irrespec- 
tive of the values of ^ , as expected from (3). For small £/’s, 


tQ / 

Z (z,to) and Z (z, to), and their corresponding power spec- 
tra F§(M,t 0 ) and F z (M,t 0 ), were almost identical. How- 
ever, as ^jf’s increased in wave fields where Kj vectors were 
not all parallel, Z (z,to) deviated more and more from 
Z (z,to) by developing nonsinusoidal waveform distortions. 
This led to significant leakage of spectral power in F§(M, to) 
to Fourier wavenumbers M > me. Eventually amplitudes 
became large enough to generate isentropic overturning in 
places, implying convective instability, turbulence and mix- 
ing. On approaching this limit, F§(M,t 0 ) attained a limit- 
ing M ~ 3 shape at M > me. 

We illustrate these general features here using a specific 
field of J = 8 waves with parameters 

mj = - ajmf , aj = 1, 1,2, 2, 3, 3,4, 4, (15) 

kj = bjkf, bj = 1,— 1,1,— 1,1,— 1,1,— 1, (16) 

C? = C/» (17) 

<i>j - 7r,7T,0,0,0,0,7r,7r, (18) 

and fundamental wavenumbers mj = 27r(10 km) -1 , kf = 
27r(50 km) 1 . From (15), the largest vertical wavenumber 
me = 4m/. An amplitude £/ =0.112 km was used for 
each wave: values only slightly larger than this produced 
regions of convective instability. Figure la plots distorted 
isentropes in the model at t = 0, using I x = 80, I z = 1000, 
Ax = 0.625 km and A z = 0.02 km. Advected parcels P ix , iz 
within a small region centered at x = 0 and z = 14 km are 
shown in Figure lb. The observational “beam” of a profiler 
at x = 0 (i x = 40) is also shown. 

Figure 2a shows Z (z, t 0 ) and z' Q ( z,t 0 ) profiles at x 4 o = 0 
and to = 0. We see a nonsinusoidal distortion of z' Q (z,to) 
compared to the sinusoidal Z (z,to) profile. Spectra (Fig- 
ure 2b) show that the distortions produce an M ~ 3 spectrum 
at M > me extending over more than a decade in M. 

The wave field (15)-(18) produces a horizontal standing 
wave pattern. At to = 0, xf 0 ,i z (to) = x 40 = 0 since x'j os- 
cillations interfere destructively here (see Figure lb). Con- 
versely, oscillations interfere constructively here, which 
yields a one-dimensional (ID) profile solution at x 40 given 
from (5), (6) and (8) by 

C (x40,Z + C (X40,Z, to),to\ = £ (X40, Z,to), (19) 



Figure 2. (a) Z (z,to) (dotted curve) and Z' Q (z, to) (solid 
curve) from model results in Figure 1 at x 40 = 0 and t 0 = 0; 
(b) Fz(Mjto) (dotted curve) and F§(M f to) (solid curve) . 
Model spectrum (23) is also plotted (thick light curve) usin g 
a = 0.4. 
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where £ (%4o,z,to) = Z (z,to). Regridding the advection 
profile £ ^X 40 ,z + £ ( X 40 , z, to), to ^ onto the regular verti- 

cal grid yields the z' Q (z,t 0 ) profile. We see from (19) 
that the nonsinusoidal “steepening” of this Z' Q (z,t 0 ) oscil- 
lation (Figure 2a) is caused by differential vertical advection 
of points on the original profile z'(z,to) (Figure lb), which 
yields a characteristic longitudinal distortion or “self advec- 
tion” of the waveform. Note from Figure lb, however, that 
both horizontal and vertical parcel displacements are im- 
portant in producing distorted Z Q (z,t) waveforms at other 
x ix ’s in the model, and for wave fields different to (15)-(18). 

Since the rrij’s are constant and have magnitudes < me, 
no identifiable waves exist in the F® (M, to) “tail spectrum” 

at M > me in Figure 2b. Indeed, on comparing Z' Q (z,to) 
and Z (z,to) in Figure 2a, we note that the wave lengths 
remain the same but the wave shapes differ. Thus, the M -3 
spectrum in Figure 2b is essentially a Fourier effect, in which 
nonsinusoidal steepening of Z (z, to) contributes a series of 
higher Fourier harmonics which asymptote to a characteris- 
tic M ~ 3 spectrum as amplitudes increase towards a “satu- 
ration” value ofC/~0.112 km. 

To illustrate further features, we utilize the ID solution 
(19) at £40, since it is conceptually useful and can be gridded 
at a high vertical resolution. Here we use I z =2 10 4 and 
Az = \ f /I z for a given m f = 27t/A/. Regridding is achieved 
here using cubic spline interpolation. A series of steepened 
wave forms similar in form to Figure 2a were created from 
(15)- (18) for arbitrary m/ by choosing 

C / — e sXf = 4e s A c, (20) 

where e s is a “saturation” constant for a given set of 0/s. 

For the phases (18), we obtained e s « 0.0115 at this 
higher vertical resolution. Left panels of Figure 3 show 



vertical displacements (km) M (eye km' 1 ) 


Figure 3. Left column: z\z,to) (dotted curve) and 
Z (z,t 0 ) (solid curve) from (20) at x 40 = 0. Right col- 
umn: spectra Fz{M , to) (dotted curve) and F^(M, to) (solid 
curve). Model spectrum (23) is also shown using a = 0.4. 
Results are for three me = m f / 4 values: (10 km) -1 (top 
row); (2.5 km) -1 (middle); (1 km) -1 (bottom). 



Figure 4. As in Figure 2, but for horizontal velocity per- 
turbations at i x = 555 (x = 2.75 km). Model “saturated” 
spectrum also plotted using a = 0.1. 


fQ t 

Z (z,to) and Z (z,to) profiles at to =0 using three m/ 
values: 2 tt( 40 km)" 1 (top row); 27r(10km)- 1 (middle row); 
m f — 27r(4km) 1 (bottom row). Corresponding spectra are 
shown to the right. F£(Af, to) spectra at M > me = 4 m f 
not only have M 3 shapes (extending over more than two 
decades in M), but their intensities are nearly identical in 
each case, despite the different m/ and values used. 

Existing spectral models (based to some extent on col- 
lated observations) yield “saturated” spectra of relative po- 
tential temperature perturbations at large M of the form 

~ a M 3 ’ ( 21 ) 

where a is a universal constant, typically in the range ~ 0.1- 
0.5 depending on the model used [Smith et al , 1987; Hines , 
1991; Gardner , 1996]. From (14), this implies a saturated 
spectrum of hydrostatic vertical displacements 

F s z a \M) = aM~\ ( 22 ) 

for M in units of rad m -1 : for units of eye m -1 , an ad- 
ditional factor (2 tt) - 2 appears on the right of (22). This 
spectrum is plotted in Figures 2b and 3 using a “best fit” 
a « 0.4 that closely matches both the shape and intensity 
of the various F^(M : to) tail spectra. Thus, even when mf 
and me vary, the advectively-generated F§ (M, to) tail spec- 
tra at M > me always converge to the “universal” M ~ 3 
spectrum observed in the atmosphere [Smith et al , 1987]. 
The “best fit” a value varied somewhat from simulation to 
simulation: e.g., adding tt/ 2 to all phases in (18) yielded 
e s ~ 0.0102 and F^(M, to) spectra at M > m c with a “best 
fit” a - 0.2. 

For a single hydrostatic wave, uj = -tN5jCj, where uj 
is the peak horizontal velocity amplitude and Sj = ±1, de- 
pending on the sign of mj/kj (assuming uj > 0) [Gossard 
and Hooke , 1975]. Combining this with (3) yields the hori- 
zontal velocity field 

j 

u {%ix > Z i z it) AT ^ ^ fijCj cos (kjXi x 4- mjZi z — ojjt -K (f>j), 

3 = 1 

( 23 ) 

which leads in the same way as before to u'(z,to) and 
U (z, to) profiles at a given x ix . The waves (15)-(18) yield 
a vanishing U (z,to) profile at x 4 o = 0, so we consider 
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profiles at other locations where both horizontal and ver- 
tical displacements will be significant (Figure lb). There- 
fore a high horizontal resolution is used here (I x = 1000, 
Ax = 50 m). 

Model-generated profiles at i x = 555 (x ix = 2.75 km) are 
shown in Figure 4a using N = 0.02 rad s -1 . Profile distor- 
tions are different here: whereas Z Q (z,to) steepened near 
the nodes of the oscillation (Figures 2a and 3), u' Q (z,to) 
departs from sinusoidality mainly around the peaks and 
troughs of the waveform. This is due to the ±90° phase 
shift between Uj and Q for each wave. Peak amplitudes of 
the profiles in Figure 4a also differ, due to nonzero hori- 
zontal parcel displacements that arise here. Spectra in Fig- 
ure 4b show formation of F§(M, t 0 ) “tail spectra” ~ M ~ 3 
at M > me • Using (22) and (23), the corresponding “sat- 
urated” horizontal velocity spectrum is « aN 2 M~ 3 [Smith 
et a/., 1987] and is plotted in Figure 4b using a nominal 
a « 0.1. Intensities of F§(M, to) tail spectra are somewhat 
less than model predictions. 


where a = (167re s ) « 0.6 for the waves (15)-(18). In an 
isothermal unsheared atmosphere, U' 2 = N 2 ^ 2 grows with 
height as exp (z/H p ), where H p is density scale height. From 
(25), me must decrease exponentially with height with a 
scale height 2H P to retain the same marginally stable pro- 
files in Figure 3. In practice, wave fields will not scale with 
height in this manner. Rather, as me decreases with height, 
those waves of |mj|~ me will be progressively removed or 
attenuated [Hines, 1993; Gardner , 1996; Eckermann , 1997], 
yielding a more gradual increase (decrease) in variance (me) 
with height [Hines, 1996]. 

/Q 

The Z (z, to) profile in Figure 2a implies large narrowly- 
confined enhancements in static stability at ~6 and 14 km. 
These features resemble the narrow “sheets” observed in 
high resolution stratospheric temperature profiles [Dalaudier 
et al. , 1994]. Turbulence often occurs near these layers, 
which also agrees with the current model, since they arise 
here only at amplitudes close to the onset of convective in- 
stability. 


Discussion 

A simple “Occam’s Razor” model reveals that the purely 
kinematic effects of a linear superposition of gravity waves 
can generate vertical wavenumber spectra similar to obser- 
vations. Simulated vertical displacement/potential temper- 
ature spectra attained limiting quasi- universal M~ 3 forms 
at large M that closely resemble observed spectra. Horizon- 
tal velocities also had M~ 3 shapes, but spectral intensities 
here were generally smaller than observed. A range of omit- 
ted effects (e.g., rotation, three-dimensionality, dynamical 
nonlinearities, refraction, different kj 1 s, spectral averaging 
etc.) may account for some of these discrepansies. 

Hines [1991, 1993, 1996] argued that advective nonlin- 
earities are important for atmospheric gravity wave spec- 
tra, citing as evidence the Allen and Joseph [1989] (AJ89) 
model of oceanic gravity waves. The AJ89 model produces 
M spectra beyond some cutoff wavenumber me- The 
present simulations, while highly idealized, have features 
which resemble the AJ89 analysis: e.g., linear superposition 
of Lagrangian wave displacements, inclusion of advection ef- 
fects, computation of both Eulerian and Lagrangian spectra, 
and M 3 “tail spectra” beyond a cutoff wavenumber me. 
Hines [1991, 1993] argues that Doppler-spread theory (DST) 
models these processes too [Hines, 1996; Chunchuzov, 1996], 
though other studies have questioned this [Eckermann, 1997; 
Broutman et al., 1997]. 

Quasi-universal M~ 3 spectra arose here due to a limiting 
“saturation” constant e s in (20) for a given set of waves. 
A similar saturation constant marks the onset of convective 
overturning for a single hydrostatic wave of wavenumber me 
and amplitude Q, 

0 = mc 1 =( 2tt) _1 Ac. (24) 

This yields Cf/\c ~ 0.16, whereas from (20) (f/\c — 

4e s « 0.05 for the simulations here (e s « 0.0115). The 

smaller latter ratio indicates that instabilities due to su- 
perposition set in first here. From (3) and (20), the total 
variance 

= i < 25 ) 
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Plate 1 

Figures 1,2,3,4,5,6,7,9,10,11,12 

(note that Figures 8,13,14 and 15 are to be found in the main body of the text) 
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In situ Meteorological Measurements Made 
on the DC-8 During the SUCCESS Campaign 
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